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ABSTRACT 

Recent discoveries of several transiting planets with clearly non-zero eccentricities and some large 
obliquities started changing the simple picture of close-in planets having circular and well-aligned 
orbits. Two major scenarios to form such close-in planets are planet migration in a disk, and planet- 
planet interactions combined with tidal dissipation. The former scenario can naturally produce a 
circular and low-obliquity orbit, while the latter implicitly assumes an initially highly eccentric and 
possibly high-obliquity orbit, which are then circularized and aligned via tidal dissipation. 
Most of these close-in planets experience orbital decay all the way to the Roche limit as the previous 
studies showed. We investigate the tidal evolution of transiting planets on eccentric orbits, and find 
that there are two characteristic evolution paths for them, depending on the relative efficiency of tidal 
dissipation inside the star and the planet. Our study shows that each of these paths may correspond to 
migration, and scattering scenarios, respectively. We further point out that the current observations 
may be consistent with the scattering scenario, where the circularization of an initially eccentric orbit 
occurs before the orbital decay primarily due to tidal dissipation in the planet, while the alignment 
of the stellar spin and orbit normal occurs on the similar timescale to the orbital decay largely due to 
dissipation in the star. We also find that even when the stellar spin-orbit misalignment is observed 
to be small at present, some systems could have had a highly misaligned orbit in the past, if their 
evolution is dominated by tidal dissipation in the star. 

Finally, we also re-examine the recent claim by Levrard et. al. that all orbital and spin parameters, 
including eccentricity and stellar obliquity, evolve on a similar timescale to orbital decay. This counter- 
intuitive result turns out to have been caused by a typo in their numerical code. Solving the correct set 
of tidal equations, we find that the eccentricity behaves as expected, with orbits usually circularizing 
rapidly compared to the orbital decay rate. 
Subject headings: planetary systems, planets and satellites: formation 



1. INTRODUCTION 

More than 450 exoplanets have been discovered so 
far. Out of about 360 extrasolar planetary systems, 
roughly 30% possess close-in planets with semimajor axis 
a < 0.1 AU. Also, there are 18 out of 45 multiple-planet 
systems with at least one close-in planet. The mean 
eccentricity for extrasolar planets with a < 0.1 AU is 
close to zero, while for planets beyond 0.1 AU, it is 
e ~ 0.25. This sharp decline in eccentricity close to 
the central star is usually explained as a result of effi- 
cient eccentricity damping due to tidal inte ractions be- 
tween the sta r and the planet (e.g., Rasio et al.l 119961 : 
iJackson et al.l [2008) . Additionally, there are currently at 
least 26 systems with measurements of the projected stel- 
lar obliquity angle A (see T able [H through the Rossiter- 
McLaughlin (RM) effect (Rossiter 1924; 'McLaughliri 
[r924 ; Ohta et al. 2005 : Gaudi fc Winn 2007) . Although 
many systems have projected stellar obliquities consis - 
tent with zero within 2a (e.g. jFabrvckv fc Winnl[2009^ . 
suggesting near-perfect spin-orbit alignment, there are 
now several planetary systems that are clearly misaligned 
IXriaud et al. 2010). Examples include XO-3, HD 80606, 
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and WASP-14, which are in prograde orbits with A ~ 
37.3 ± 3.7, 53I211 a-nd — 33.1 ± 7.4 degrees, respectively 
(jWinn et al.l l2009cl lal: iJohnson et all l2009( l. as weU as 
HAT-P-7, WASP-2, WASP-8, WASP-15, and WASP- 
17, which have retrograde orbits with A ~ 182.5 ± 9.4, 
-1531:11, -1 20±4, -139.61^;^ and -147.3J: p degrees, 
respectively (jWinn et al.ll2009bl : PTviaud et al.l l2010). 

The standard planet formation theory predicts that 
giant planets are formed beyond the so-called "ice line" 
(a > 3 AU) , where solid material is abundant due to 
condensation of ice. To explain the orbital properties of 
close-in planets, two different scenarios have been pro- 
posed. Both of them could potentially explain the prox- 
imity of these planets to the stars, but they predict dif- 
ferent distributions for orbital inclinations, and possibly 
eccentricities. One scenario is orbital migration of plan- 
ets due to gravitational interactions with gas or plan- 
etesimal disks (e.g. , iGoldreich fc Tremaind Il980l : IWardI 
Il997t iMurrav et al.lll998[ ). which would naturally bring 
planets inward from their formation sites. The scenario 
can also account for the observed small eccentricity and 
obliquity seen in the majority of close-in planets, be- 
cause the disks tend to damp e ccentricity and inclina- 
tion of the planetary orbit (e.g.. | Goldreich fc Tremaind 
119801 : iPapaloizou fc Larwoodll2000( ). Alternativelv. such 
close-in planets can be formed by tidally circularizing 
a highly eccentric orbit. A natural way of initially 
increasing the orbital eccentricity is through gravita- 
tional interactions between several planets. Although 
such interactions alone may not be able to populate 



the inner re gion of planetary system s (less than 0.1 — 
1 AU. e.g.. lAdams fc LaughlinI l2003t IChatteriee et all 
l2008t iMatsumura et ahl |201CID . the orbital eccentricity 
can be increased due to scatter ing, ejection, or Kozai 
cycles (jKozail Il962t iLidovl [196^ , so that the pericen- 
ter of the planetary orbit becomes small enough for 
tidal inter actions with the central star to become impor- 
tant fe.g., iNagasawa et al.l 120081 ). These gravitational 
inter actions also tend to in crease the orbital inclina- 
tion. IChatteriee et al.l (|2008^ performed a number of dy- 
namical simulations of three-planet systems, and showed 
that the final mean inclination of planetary orbits is 
about 20 degrees, and that some planets could end up 
on retrograde orbits (about 2%, S. Chatterjee, private 
communication). When Kozai cycles increase the or- 
bital eccentricity, the process is called Kozai migration 
l|Wu fc Murravll2003| ). Kozai migration occurring in bi- 
nary systems ma y be responsible for at least some of the 
close-in planets (iFabrvckv fc Tremaind I2007L IWu et all 
I2007t iTriaud et al.ll2010 f). One of the goals of our study 
is to explore the possibility and implications of forming 
close-in planets via tidal circularization of a highly ec- 
centric planet. 

Independent of their formation mechanism, these close- 
in planets are currently subject to strong tidal interac- 
tions with the central star, and such interactions could 
dominate the orbital evolution of these planets. In multi- 
planet systems, secular planet-planet interactions may 
also affect the orbital evolution. However, in Section 2, 
we will show that it is unlikely that the current and fu- 
ture evolution of the observed close-in planets is strongly 
affected by any known or yet-to- be-detected companion 
(see also iMatsumura et al.ll2008[) . In this paper, we in- 
vestigate tidal evolution of close-in planets to distinguish 
the two formation scenarios of them. 

Tidal evolution in a two-body system leads to either a 
stable equilibrium state, or to orbital de cay all the wa y 
to the Roche limit ("Darwin instability" iDarwinI I1879D . 
Such a stud y for e xoplanetary systems was first done by 
iRasio et"aLl (|1996[ ). who suggested that 51 Peg b, the 
only close-in planet known at the time, would b e Dar- 
win unstable. In a recent paper, iLevrard et al.l (|2009l 
hereafter, LWC09) investigated the tidal evolution of all 
transiting planets, and pointed out that most of these 
planets (except HAT-P-2b) are indeed Darwin unstable, 
and thus undergo continual orbital decay, rather than ar- 
rive at a stable, equilibrium orbit (see Section 4.1, and 
Table [2] for updated results). These Darwin- unstable 
planets may eventually be accreted by the central star, 
which ha s been suggested fo r some systems obs ervation- 
ally (e.g. [Gonzalez 1997; E cuyillon et al.ll2006[ ) and nu- 
merically (I Jackson et al.ll2009l) . 

There is some confusion in the literature concerning 
the evolution timescales for the various orbital elements. 
LWC09 studied tidal evolution of transiting planets by 
taking into account energy dissipation in both the star 
and the planet. They concluded that all orbital and spin 
parameters, with the exception of the planetary spin, 
would evolve on a timescale comparable to the orbital 
decay timescale. This would imply that both circularity 
of the orbits and spin-orbit alignment seen in many sys- 
tems are primordial, because neither obliquity nor eccen- 
tricity could be damped to zero before complete spiral- 
in and destruction of the planet. Unfortunately, after 



much investigation and comparison with their work, we 
determined that LWC09 had a typo in their code (con- 
firmed by B. Levrard, personal communication), which 
made them underestimate the energy dissipation inside 
the planet by several orders of magnitude. By integrating 
the correct set of tidal evolution equations, we find that 
there are two characteristic evolutionary paths depend- 
ing on the relative efficiency of tidal dissipation inside the 
star and the planet. When the dissipation in the planet 
dominates, the eccentricity damping time is shorter than 
the orbital decay time {t^, ^ r^ ; see Se ction 4 for details) , 
exactly as expected intu i tively (e.g., iRasio et al. Ill 9961 : 
iDobbs-Dixon et al.l[200llMardling fc Linll2004l ). On the 
other hand, when the dissipation in the star dominates, 
the eccentricity damps on a similar timescale to the or- 
bital decay. We will show that the latter path is funda- 
mentally different from what is suggested by LWC09 in 
Section 4.2. 

There have been many other recent s tudies of tidal 
evolution for exoplanets. 'Jacks on et al.l (|2008") empha- 
sized the importance of solving the coupled evolution 
equations for eccentricity and semi-major axis. Integrat- 
ing their tidal evolution equations backwards in time, 
they showed that the "initial" eccentricity distribution of 
close-in planets matches more closely that of planets on 
wider orbits; they suggested that gas disk migration is 
therefore not responsible for all the close-in planets. We 
reconsider th e validity and imp l icatio n of such a study 
in Section 5. iBarker fc Ogilvid (|2009( l studied the evo- 
lution of close-in planets on inclined orb its, by includ- 
ing t he effect of magnetic braking (Dob bs-Dixon et all 
l2004f) . and pointed out that a true tidal equilibrium 
state is never reached in reality, since the total angu- 
lar momentum is not conserved due to magnetized stel- 
lar winds. They also showed that neglecting this effect 
could result in a very different predicted evolution for 
the systems they considered. Throughout this paper, 
we compare the tidal evolutions with and without the 
effect of magnetic braking. Another potentially impor- 
tant effect caused by tidal dissipation is the i nflation of 
the p lanetary radius (Bodenhcimer ct al. 2001: JGu et al.l 
|2003( ). Many groups have explored this possibility, mo- 
tivated by observa tions of infla ted radii for some tran- 
siting planets (e.g . iBarge et al.l[20 08: Aloi iso et alll2008l : 
iJohns-Krull et al.l 120081 : iGillonet al.i i200i). It appears 
that at least some of these inflated planet s could be ex- 
plained as a result of past tidal heat ing (Jac kson et al.l 
[2008t IM ller"dran[2009|: llbgui fc Burrows .2001) . More 



recently. iLeconte et al.l (j2010f ) revisited this problem, and 



pointed out that the truncated tidal equations used in 
many previous studies could lead to an erroneous tidal 
evolution for moderate to high eccentricity (e > 0.2). 
Solving the complete set of tidal equations, they showed 
that the orbital circularization occurs much earlier than 
previously estimated, and thus only moderately bloated 
hot Jupiters could be explained as a result of tidal heat- 
ing. In this paper we neglect this effect entirely and treat 
the planetary radius as constant for simplicity (but see 
Section 5.1). 

The outline of the paper is as follows. First, we jus- 
tify our approach in Section 2 by showing that the cur- 
rent/future evolution of these planets is likely dominated 
by tidal dissipation, rather than by their gravitational in- 
teraction with a more distant object (planet on a wider 



orbit or distant binary companion). Then we present 
our set of tidal evolution equations in Section 3. We 
re-examine the tidal stability of transiting planets, and 
investigate the tidal evolution forward in time in Sec- 
tion 4. We identify two characteristic evolution paths for 
Darwin-unstable planets, and also discuss the results of 
LWC09. In Section 5 we explore the past evolution of 
transiting planets and its implication. Our results im- 
ply that each evolutionary path may be consistent with 
migration, and gravitational-interaction-induced forma- 
tion scenarios of transiting planets, respectively. We also 
point out that in a limited case, it's possible to have a sig- 
nificant stellar obliquity damping before the substantial 
orbital decay to the Roche limit. In Section 6, we study 
the different definitions of tidal quality factors, and in- 
vestigate their effects on evolution. Finally, in Section 7, 
we discuss and summarize our results. 

2. POSSIBLE COMPANIONS TO CLOSE-IN PLANETS 

In this section, we assess the importance of a 
known/unknown companion on the current and future 
evolution of close-in transiting planets. Secular or res- 
onant interactions with other planets or stellar com- 
panions could potentially perturb the planetary orbits 
significantly. For example, when there is a large mu- 
tual inclination between their orbits (> 39.2°), Kozai- 
type (quadrupole) pertur bations can become important 
(|Kozailll962t lLidovlll962f ). Such highly misaligned sys- 
tems may naturally occur for binary syst ems with semi- 
major axis > 30 - 40 AU JHald 119941), or as a re- 
sult of planet-planet scattering (jChatteriee et all 120081 : 
iNagasawa et al.ll2008[ ). For smaller mutual inclinations, 
octupole-level perturbations may still moderately excite 
the orbital eccentricity, as long as the companion has a 
non-circular orbit. 

For this secular perturbation from a companion to be 
affecting the current and future evolution of close-in plan- 
ets, it must occur fast compared to other perturbations 
that cause orbital precession. These competing per- 
turbations include general relativistic (G R) precession, 
tides, as weh as rotational distortions (jHolman et al.l 
Il997t ISternellT939D . Since the effects of pericenter pre- 
cession due to stellar and planetary oblateness are usu- 
ally small comp ared to those caused by GR prece ssion 
iKiselevaet all 119981: iFabrvckv fclVemainel I2007D . we 
simply neglect rotational distortions here. 

The pericenter precession timescales corresponding to 
GR, Kozai-type perturbations (due to a high-inclination 
perturber), and secular coupling to a low-inclination 
pertu r ber, can be written as fol l ows (jKiseleva et aLl 



c is the speed of light. In r, 



19981 IFabrvckv fc Tremah^ [20071 iZhou fc SunI 120031 : 
Takeda et al.ll2008D : 
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with the standard Laplace coefficients h\,r^{a-p/ac) {i = 
1,2). For the secular timescale, the upper sign is cho- 
sen when Mp < Mc and the lower one is chosen when 
Mp > Mc- Such an appr oximation is reason able when 
the system is hierarchical (jTakeda et al.|[2008l ) and thus 
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Following the approach of iMatsumura et al.l ([2008') , we 
compare these various precession timescales to determine 
which perturbation dominates. 

Fig. [1] shows the resulting constraints on the mass and 
orbital radius of the hypothetical planetary/stellar com- 
panion for CoRoT-7b and HAT-P-13b. These are the 
only two known multi-planet systems with a transiting 
planet, and the companions are shown as the filled circles 
in the plot. In order to perturb the inner planet signif- 
icantly despite the GR precession, the companion must 
exist left of the blue lines to induce Kozai oscillations, 
and left of the orange line for octupole perturbations to 
be important. The kink seen in the orange line occurs 
where Mp '^ Mc and thus our approximation in calculat- 
ing the secular timescale breaks down. For HAT-P-13, 
we find that HAT-P-13c is located left of the blue lines, 
but right of the orange one. Thus, we expect that HAT- 
P-13c can perturb the orbit of HAT-P-13b significantly 
only if they have a large mutual inclination of > 39.2°. 
For CoRoT-7, on the other hand, we find that the secu- 
lar timescale due to CoRoT-7c is compa rable to the GR 
timescal e (r-np ~ tor). As show n by Fo rd et all (j2000l ) 
(see also I Adams fc LaughlinI 2006), such a resonance can 
increase the eccentricity significantly. However, for small 
eccentricities e ~ 0.001, we do not find any significant ec- 
centricity oscillation below a mutual inclination of about 
40°. Thus it is possible that the eccentricity of CoRoT- 
7b is significantly oscillating if the misalignment of the 
orbit of CoRoT-7c with respect to CoRoT-7b is large. 

In Fig. [51 we present similar plots for all transiting 
planets on an eccentric orbit. Vertical and horizontal 
lines are drawn for reference at 1 AU and lOAf j, respec- 
tively. It is clear that all systems except HD 80606 are 
unlikely to have a companion which has a secular per- 
turbation timescale shorter than the GR timescale (left 
of orange and blue lines), and at the same time does not 
cause detectable radial-velocity variations (below black 
dotted lines, or beyond the observation limit in semi- 
major axis). The figure predicts that, if a hypothetical 
planet which can cause a significant secular perturbation 
is too small to be observed (i.e., a planet exists left of 
the orange line and below black dotted lines), its mass 
would be comparable to Earth or smaller. The figure also 
predicts that, if such a planet cannot be observed now 
simply because we are not observing long enough (i.e., a 
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Fig. 1. — Mass and semi-major axis of a hypothetical compan- 
ion which can affect the orbital evolution of the observed close-in 
planets in CoRoT-7 and HAT-P-13. Orange line indicates the re- 
gion within which the secular timescale becomes short compared 
to the GR timescale. The kink occurs where the approximation 
for the secular timescale breaks down (see text). Left and right 
blue dashed lines indicate similar boundaries for Kozai cycles with 
the companion eccentricity of 0.2, and 0.9, respectively. Black dot- 
ted lines indicate the radial velocity detection limits for 3, 30, and 
100 m/s. Both of these systems have a known second planet, which 
is indicated by a filled circle. 
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Fig. 2. — Similar plots to Fig.[T]for all transiting planets with an 
eccentric orbit. Vertical and horizontal lines are drawn at 1 AU and 
lOA/j for comparison. It is clear that all systems except HD 80606 
are unlikely to have a companion which can cause significant secular 
perturbations (faster than the GR precession). 



planet exists left of the orange line, and beyond, for ex- 
ample, 1 AU) , its mass would be comparable to brown 
dwarfs or larger. 

HD 80606b is in a wid e (- 1200 AU) binary system 
(jEggenberger et al]|2004D . However, as can be seen in 
the figure, this companion cannot induce secular pertur- 
bations fast enough compared to the GR precession. 

Resonant interactions can work in a similar way, but 
it is unlikely that all of these planets have such a com- 
panion. Thus, the evolution of currently observed tran- 
siting planets with an eccentric orbit is likely dominated 
by tidal dissipation rather than interactions with outer 
companions. 

3. TIDAL EVOLUTION EQUATIONS 

In this paper, we numerically study the evolution of 
observed transiting planets by integrating a set of equa- 
tions describing the tidal interactions, and by assuming 



that the known/unknown companions' effects can be ne- 
glected. 

We follow the general approach of the equilibri um tide 
mode l with the weak friction approximation ([Darwinl 
Il879( ). The effects of the dynamical tide are neglected 
for simplicity. In absence of any dissipation, the tidally 
distorted body is assumed to take the equilibrium shape 
that adjusts itself to the external potential field of the 
tide-raising body. When the dissipation is non-negligible, 
the equilibrium surface is either lagged or led, depend- 
ing on whether the spin frequency is smaller or larger 
than the orbital frequency. In the limit of small viscosi- 
ties, this phase lag {(j)) can be approximated to be pro- 
portional to the tidal forcing frequency {a) as ^ Aicr. 
This allows us to interpret the phase lag as the tidal bulge 
that could have been raised a constant time Ai ago in an 
inviscid case. Within the context of this model, we can 
derive the secular equations which are valid for any value 
of eccentricity and o bliqu i ty, by following the approach 
of lAlexanded (|1973[ )-. iHutI (|198li ) (also see iLeconte et all 
120101 ). Taking into account tides raised both on the cen- 
tral star by the planet and on the planet by the star, the 
complete set of tidal equations can be written as follows 
for the semi-major axis a, eccentricity e, stellar obliq- 
uity e*, planetary spin lo^, and stellar spin w^,, respec- 
tively (e.g.. Hut 1981; Levrard e t al. 2007). Note that 
we assume that the planetary obliquity is zero (i.e., the 
equatorial plane of the planet coincides with the orbital 
plane). We explicitly write down all equations in order 
to compare our results with UWC09: 
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The subscripts * and p denote the star and pl anet, 
resp ectively. These equations agree with those of iHutI 
ifT^lt ) in the limit of small e*. In LWC09, -q = 

"* ^Mp^'' $(l - e^)"^^^V is set to zero (B. Levrard, 
private communication). This term 77 is smaller than 
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COS e* for most systems, but can be comparable to, or 
larger than cose* for some systems. Examples include 
CoRoT-1, HD 149026, Kepler-4, Kepler-8, and WASP- 
17. In the above equations, k2 is th e Love num ber for 
the second-order harmonic potential (JLovgllMJ), Ai is 
a constant time lag, n is the mean motion, and a is the 
square of the radius of gyration with a* = 0.06, and 
ap = 0.26. The eccentr icity f unctio ns /i(e^) — /5(e^) are 
defined as follows as in iHutI (|198lD . 
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Although it is important to solve the coupled 
semima.jor axis and eccentricity evolution equations 
(| Jackson et al.l[2008[ ). the semi-major axis evolution of 
a planetary system is likely dominated by the energy 
dissipation in the star. This is because, in an eccen- 
tric orbit, a gaseous planet's rotation will be tidally 
damped to an asymptotic state that is somewhat faster 
than a value synchronous with the orbital mean motion. 
The tidal torque is strongest at pericenter where the or- 
bital angular velocity exceeds the orbital mean motion, 
and therefore the tidal torque averaged around the or- 
bit vanishes when the rotation rate exceeds the mean 
motion n. This asymptotic state is often referred to as 
pseudo-synchronous rotation. When the planetary spin 
approaches pseudo-synchronization with the orbit nor- 
mal ujp ^ n, the contribution from the second term in 
Eq. [7] becomes negligible for planets with small eccen- 
tricities. Therefore, unless the eccentricity is very high, 
the orbital evolution is largely determined by tidal dissi- 
pation in the star. 

It's not immediately clear whether the tidal energy dis- 
sipation leads to either orbital decay or orbital expansion. 
For all transiting systems except CoRoT-3, CoRoT-6, 
HAT-P-2, HD 80606, and WASP-7, the host star rotates 
slowly compared to the orbit (i.e., w, < n, see Table [5]). 
Therefore, tidal dissipation in the star tends to lead to or- 
bital decay, by transferring angular momentum from the 
orbit to the stellar spin. On the other hand, when the 
host star is rapidly spinning, planets could migrate out- 
ward, which may have prolonged the lifetime of some of 
the exoplanets (Dobbs- Dixon et al.ll2004( ). In Section 4.1, 
we show that HAT-P-2 is Darwin unstable and migrating 
inward, while the other systems with a rapidly rotating 
host star are evolving toward the stable tidal equilib- 
ria. More specifically, CoRoT-3, CoRoT-6, and WASP-7 
are currently migrating outward toward the stable tidal 
equilibria, while HD 80606 is migrating inward toward 
the stable state. 

By comparing Eq. [7|and[51 we see that the eccentricity 
may be damped on a similar timescale to the semima- 
jor axis, when the eccentricity damping is dominated by 
tidal dissipation in the star (i.e., the first term in Eq. [5] 
is much larger than the second term). We discuss this 



further in the next section. 

In some of our calculations we also take into account 
the stellar magnetic braking effects, and thus the loss 
of angular momentum due to stellar winds. For this 
purpose, we assume that the Skumanich's Law describes 
the decrease of the stellar spin sufficiently well so that 
the average surface rotation velocities of stars that are 
not interacting with close-in planets can b e related to 
the stellar age (rage) as F, sini* oc 1/ ^ff^ (jSkumanichl 
Il972l ). From this relation, the change in stellar spin can 
be written as follows. 
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where 7 is a calibration factor, and the subscript de- 
notes the normalization factors. By choosi ng Kto = 
4 km/s and Tagc, = 1 Gyr as in Dobbs-D ixon et al.l 
(|2004f l. we can define {3 = 7I.5 x lO^^^^yr. We adopt 
7 = 0.1 for F stars, and 7 = 1 for G, K, and M stars 
as in [Barker &: Qgilvii (|2009f ). For F stars, the smaller 
calibration factor is chosen since magnetic braking is less 
efficient due to the very thin, or completely absent outer 
convective layer. We add this —fiuj^ to Eq. [TT] when 
including the effects of magnetic braking. 

3.1. Tidal Quality Factors 

It is common to describe the dissipation inefficiency 
of tides in terms of the tidal quality factor instead of a 
constant lag angle 0, or a constant time lag Ai. The spe- 
cific dissipation function is defined as follows (|Goldreichl 
Il96l 

2itE* I 

§- (dE/dt) dt ^ t^ ^ ' 

where E* in the numerator is the peak energy stored 
in tides during one tidal cycle, while the denominator 
represents the energy dissipated over the cycle. When 
the phase angle (j> (which is twice the geometrical lag 
angle) is small, this is simplified as Q ~ l/i^, which 
implies a large Q value. The estimated tidal qual- 
ity factors are > 10'^ for Jupiter a nd Neptune (e.g., 
iLainev et al.l 120091 ; iZhang fc Hamilto n 2008) , and even 
larger val ues are expected for synchronized close-in ex- 
oplanets ([Qgilvie fc Lin|[2004D . On the other hand, the 
usua lly adopted va.lues a re > 10^ for main sequence stars 
(e.g. iTrilling et al]|1998D . Thus, we are generally inter- 
ested in the case of Q ^ 1, and the above approxima- 
tion is reasonable. Using the weak friction approximation 
(0 ~ Ata) , we can redefine Q as 



Q 



TAt 



(14) 



In the following sections, we study some limiting cases 
of tidal frequencies. Before the spin-orbit synchroniza- 
tion, the tidal dissipation is generally dominated by the 
semi-diurnal tide with the forcing frequency of ct = 
\2uj-2n\ (Ferraz-Mello et al. 2008). In this case, the 
tidal quality factor can be written as 



Q 
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At\2u}- 2n\ 



(15) 



As the system approaches synchronization, the effect due 
to the semi-diurnal tide diminishes, and the annual tide 
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with a = \lijj-n\ prevails (|Ferraz-Mello et all [2008D . 
The corresponding tidal quality factor can be written 
in a similar manner to the above as 
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Ai|2w- 
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Note that when uj = n, the most efficient energy dissipa- 
tion occurs on the same timescale as the orbital period. 

In Section 4 and 5, we rewrite eg . [TllTllbv assuming that 
both stellar and planetary tidal quality factors change as 
Q DC 1/n, while in Section 6, we adopt Q* ex 1/ \2uj — 2n\ 
and Qp oc 1/ |2aj — n|. The former is chosen to compare 
our results with the study of LWC09. For the latter, 
we are implicitly assuming that close-in exoplanets are 
pseudo-synchronized, while their host stars are not. We 
will see that this is a reasonable approximation in Sec- 
tion 4. 

Throughout this paper, we discuss our results in terms 
of the modified tidal quality factor Q' = 1.5Q/fc2, where 
k2 is the Love numb er for the second-order harmonic 
potential (|Lovelll944l ). 



4.1. 



FUTURE TIDAL EVOLUTION OF TRANSITING 
PLANETS 

Two Evolutionary Paths for Darwin- Unstable 
Extrasolar Planets 



In this subsection, we revisit the tidal stability prob- 
lem for transiting extrasolar planets, and we show that 
there are two distinctive evolutionary paths for Darwin- 
unstable systems. Throughout this subsection, we as- 
sume that the total angular momentum is conserved (i.e., 
magnetic braking is neglected). We take into account the 
magnetic braking effects later in Section 4.3. 

The existence and stability of tidal e quilibrium states 
were i nvestigated by many authors fe.g.[Hutlll980l : iPeald 
Il986t IChandrasekha^ Il987t iLai et all Il99l for binary 
systems and the Solar System. Minimizing the total en- 
ergy under t he constrai nt of conservation of total angular 
momentum, iHutI (|1980t ) found that all equilibrium states 
are characterized by orbital circularity, spin-orbit align- 
ment, and synchronization of the stellar rotation with 
the orbit. He also showed that equilibrium states ex- 
ist only when the total angular momentum of the system 
Ltot is larger than some critical value Lcrit, and that such 
equilibrium states are unstable when the orbital angular 
momentum is less than 3 times the total spin angular 
momentum (ispin/-^orb > 1/3)- In this paper, we are 
interested in the ultimate fate of close-in exoplanets, and 
thus we call a system "Darwin-stable", only when it is 
expected to evolve ultimately toward a stable equilibrium 
state. All other systems are called "Darwin-unstable". 

First, we check whether the stable tidal equilibrium 
states exist for the currently known transiting planets 
listed in Table [T] Tidal equilibrium can only exist when 
both primary and secondary, with zero obliquities, are 
synchronously rotating with their orbital motion. Under 
these constraints, the total angular momentum 
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Here, C = aME? is the moment of inertia, and Lorb = 
, * '' \/Ga(l — e^) is the orbital angular momentum. 

For Ltot < icrit there can be no tidal equilibrium; for 
Ltot > icrit two equilibrium states exist. The inner state 
(Loi-b < 3Lspin) is unstable, so the only stable tidal equi- 
librium's theoutei_statejWhich requires Lorb > 3isp,;„ 
(e.g.. iHutI [19801 IPeald 119861 also see the top panel of 
Fig. [3]). A local example of du al synchrono us rotation is 
the Pluto-Charon system (e.g.. lPealg|1986[) . 

The results are summarized in Table [51 As already 
pointed out by LWC09, most systems have no tidal equi- 
librium states (i.e., Ltot/^crit < 1), and thus are Darwin 
unstable (i.e., the planet will eventually fall all the way 
to the Roche limit of the central star, even when the to- 
tal angular momentum is strictly conserved ) . Not e that, 
using the refined parameters in iPal et all (|2010( ). even 
HAT-P-2, which was the only system with tidal equilib- 
ria in LWC09, in fact now appears to have no such states 
(itot/icrit = 0.995). 

There are several systems which have tidal equilib- 
ria (itot/^crit ^ 1), and thus could be Darwin-stable, 
namely, CoRoT-3, CoRoT-6, HAT-P-9, HD 80606, 
Kepler-8, WASP-7, and WASP-17. Among these, 
CoRoT-3, CoRoT-6, and HD 80606 are clearly evolv- 
ing toward the stable equilibrium state since they all 
have L*.spin/^orb < 1/3. Therefore, these systems are 
Darwin-stable. Out of the systems with -Z^^^spjn/Lorb > 
1/3, WASP-17 is Darwin-unstable independent of the 1/3 
criterion, since it is in a retrograde orbit (Lorb < 0) and 
thus has Lgpin > Lorb- For the other systems (HAT-P- 
9, Kepler-8, and WASP-7), we checked their tidal equi- 
librium states as in Fig. [3] Here, we compare the to- 
tal angular momentum and total energy for dual syn- 
chronous state with the current values. As clearly seen 
in the figure, HAT-P-9 is Darwin-unstable, since the sys- 
tem is inside the inner unstable equilibrium state. Most 
of the angular momentum in the system is in the spin, 
and the planet falls toward the central star as a result of 
tidal energy dissipation. Similarly, Kepler-8 is Darwin- 
unstable since the system is far inside the inner unstable 
equilibrium state. On the other hand, WASP-7 exists 
just outside the inner unstable equilibrium state, and 
thus is migrating outward, toward the stable tidal equi- 
librium state. Similar plots for CoRoT-3, CoRoT-6, and 
HD 80606 reveal that CoRoT-3 and CoRoT-6 are migrat- 
ing outward toward the stable equilibria, while HD 80606 
is migrating inward toward the stable state. For systems 
with non-zero eccentricities (CoRoT-3, HD 80606, and 
WASP-17), we show the evolution explicitly in Fig. [S| 
In short, we find that all planetary systems in Table [T] 
except CoRoT-3, CoRoT-6, HD 80606, and WASP-7 are 
Darwin-unstable. 

The rest of the transiting systems have no tidal equi- 
librium (Ltot/Lcrit < 1) with the fiducial orbital and 
spin parameters listed in Table [TJ and thus should be 
Darwin-unstable. However, some of them are borderline 
systems with Ltot/^crit — 1; and they could be either 
Darwin-stable or -unstable within the observational un- 
certainties or depending on the value of the unknown 
parameters (e.g., stellar obliquity). As an example, we 
show the tidal evolution of a hypothetical planet with 
Mp — 3Mj and Rp = 1.2Rj orbiting a Sun-like star. We 
assume the initial semimajor axis, eccentricity, and stel- 
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Fig. 3. — Tidal equilibrium curves as a function of orbital 
separation for HAT-P-9. The total angular momentum and to- 
tal energy for dual synchronous states are plotted in blue curves, 
while the corresponding curves with the constant angular momen- 
tum are plotted in orange. The current values are indicated as 
the star, while open and filled circles with u and s denote unsta- 
ble, and stable tidal equilibria, respectively. Although HAT-P-9 
has tidal equilibria, the system currently exists inside the unstable 
tidal equilibrium state. Therefore, HAT-P-9 is Darwin unstable 
and migrates toward the central star as the energy dissipates. 



lar velocity oi a — 0.06 AU, e = 0.3, and w* — lOkm/s. 
Such a system would have a stable tidal equilibrium if 
the stellar obliquity is small, because Ltot/-^crit ~ 1.040 
and i»,spin/-^orb ~ 0.193 < 1/3. In Fig. [j] we show the 
results of tidal evolution for two different stellar obliqui- 
ties of e^: = 20° and 60°. Here the tidal quality factors 
are scaled as Q' = QqUq/u, where indicates the ini- 
tial/current values. This scaling ensures that our model 
is consistent with the constant time lag model as shown 
in Section 3.1. For each obliquity, we performed three 



different runs with the same Q'^ g = 10^ and different 
Q'po of 10^ 10^, and 10^. For the smaller stellar obliq- 
uity (e* = 20°), we find that the system arrives at the 
stable tidal equilibrium as circularization, synchroniza- 
tion, and alignment are achieved. On the other hand, 
for the larger stellar obliquity (e, = 60°), the system 
turns out to be Darwin unstable, and the planet spirals 
into the star on a ^ 10 Gyr timescale. The examples 
of these borderline systems include HAT-P-2, WASP-10, 
and XO-3. XO-3 was a Darwin-stable system with previ- 
ously obtained observed parameters, while HAT-P-2 and 
WASP-10 can be Darwin-stable within the uncertainties 
as we see in Fig. [B] 

The tidal evolution is not completely simple even for 
clearly Darwin-unstable systems. The bottom panels of 
Fig. m demonstrate that such systems can take either of 
two different evolutionary paths, depending on the rel- 
ative efficiency of tidal dissipation inside the star and 
inside the planet. With the smaller planetary tidal qual- 
ity factors and thus with more efficient tidal dissipation 
in the planet {Q'q = 10^, and 10^ in the figure, which 
correspond to dotted and solid curves, respectively), the 
planetary orbit circularizes before the planet spirals into 



the Roche limit (rg < Tq), while with the larger plane- 
tary tidal quality factors {Q' q — 10^ in the figure, cor- 
responding to dashed curves), the circularization time 
becomes comparable to the orbital decay time (re ~ Ta). 
The difference occurs because the orbital decay time is 
largely determined by dissipation inside the star, while 
the circularization time can be determined by either dis- 
sipation inside the star, or that inside the planet. This 
is apparent from the figure — the orbital decay times 
are similar for different Q' q values, which means that 
the semi-major axis evolution is largely independent of 
tidal dissipation in the planet, and instead is determined 
entirely by dissipation in the star (r^ ~ Ta,*)- On the 
other hand, the eccentricity damping timescales change 
significantly for different Q' q values, which suggests that 
dissipation in the planet plays a significant role in cir- 
cularization. However, the circularization time is never 
longer than the orbital decay time. Therefore, when the 
expected circularization time from Q' is longer than the 
orbital decay time, the eccentricity damping time is also 
determined by the stellar tidal dissipation (re ~ Te_*, 
which corresponds to Q' q = 10^ case in the figure). We 
show that these results are generally true in Section 4.3. 

Unfortunately, it is nontrivial to constrain tidal qual- 
ity factors and determine which type of evolution each 
system would follow. This is because tidal quality fac- 
tors depend sensitively on the detailed interior structure 
of the body (either star or planet), as well as the tidal 
forcing frequency and amplitude, and are unlikely to be 
expressed a s a simple c onstant value (e.g.. lOgilvie fc Lig 
12004 I2007t IWul |2005[ ). However, since we observe a 
sharp eccentricity decline within a < 0.1 AU, and many 
Darwin-unstable extrasolar planets are observed to be 
on nearly circular orbits, it is likely that the eccentricity 
damping time tends to be short compared to the orbital 
decay time (te < Ta). We discuss this issue further in 
Section 5.2. 

In summary, we find that there are two evolution- 
ary paths for Darwin- unstable planets — the "stellar- 
dissipation dominated" case in which all the parame- 
ters except planetary spin evolve on a similar timescale 
(te ~ Tq ^ Tf^ ^ Tu._,), and the "planetary-dissipation 
dominated" case in which circularization occurs before 
any significant orbital decay (re < Ta ''-^ t^, ~ ''"cj. )■ In 
the former case, not only the evolution of stellar spin and 
obliquity, but also that of semi-major axis and eccentric- 
ity are controlled by the efficiency of tidal dissipation 
in the star. In the latter case, the eccentricity damping 
is driven by dissipation in the planet, while the orbital 
decay is largely determined by dissipation in the star. 



4.2. Comparison with lLevrard et ali l200A ) 

Note that, although the stellar-dissipation-driven case 
(te ~ Ta) looks similar to the one illustrated by LWC09, 
there are fundamental differences. In Fig. 2 of UWC09, 
the exponential eccentricity damping approximation, in 
which they only integrate the planetary dissipation term 
in the eccentricity evolution, shows a faster damping 
time than the tidal evolution involving both stellar and 
planetary energy dissipation. Generally, however, this 
approximation should provide a timescale comparable 
to, or longer than, what is obtained with the full in- 
tegration. This is because most planet-hosting stars are 
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Fig. 4. — Tidal evolution of a planet with 3Mj at 0.06 AU and 
e = 0.3. In the left panels, black, and orange curves show the 
evolution of semi-major axis and eccentricity, respectively. In the 
right panels, orange curves show the evolution of orbital period, 
while blue and red curves show that of planetary and stellar spin 
periods, respectively. The black curves indicate the evolution of 
stellar obliquity. Dotted, solid, and dashed curves correspond to 
cases with the same Q'^ g = 10^ and different Q' q of 10^, 10®, 

and 10 , respectively. Top panels show the cases with e* = 20°, 
which are Darwin-stable, and bottom panels show the cases with 
e, = 60° , which are Darwin-unstable. 



spinning slower than the orbital frequency (see Table [5]), 
and thus energy dissipation in the star accelerates the 
ecce ntricity damping, rather than slowing it down (see 
also iDobbs-Dixon et al.ll2004l ). Therefore, for most sys- 
tems, the exponential damping approximation involving 
only the dissipation in the planet should provide an upper 
limit for eccentricity damping. 

We demonstrate this in Fig. \5\ by taking HD 209458 
(the same system considered in LWC09) as an exam- 
ple. Note however that, with current data, HD 209458 
has an orbital eccentricity consistent with zero, and thus 
is not included in the analysis presented in the rest of 
this paper. We adopt the same initial conditions and 
assumptions, and use the same set of tidal equations 
as LWC09. Here the tidal quality factors are initially 
Qp,o = Q'*,o — 10^, and scale as Q' — QqUo/u. Our re- 
sult is shown in the top panel of Fig. \E\ We find that the 
tidal evolution in Fig. 2 of LWC09 is recovered except for 
the eccentricity. In our results, the eccentricity evolves 
on a timescale consistent with the exponential eccentric- 
ity damping approximation (dashed line, which is com- 
pletely overlapped with solid line in the top left panel). 
In other words, the eccentricity damps on the timescale 
determined by tidal dissipation in the planet, as previ- 
ously shown by many authors (e.g., Rasio ct al. 1996; 
iDobbs-Dixon et al.|[200llMardling fc Linll2004D . For the 
middle panels, we assume a less efficient tidal dissipation 
for the planet {Q'q — 10^). In this case, we obtain an 
eccentricity evolution similar to LWC09, but the plane- 
tary spin-orbit synchronization, as expected, occurs more 
slowly. Note that the dashed curve is the eccentricity 
damping approximation with Q' q = 10^. Finally, for 
the bottom panels, we reproduce Fig. 2 of LWC09 by 
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Fig. 5. — Evolution of the orbital and spin parameters of 
HD 209458 (cf. Fig. 2 of LWC09). Three different cases are shown 
from top to bottom. In all cases, the integrations are stopped when 
the planet formally hits the stellar surface. In the left panels, black 
and orange curves represent the evolution of semi-major axis and 
eccentricity, respectively. The dashed orange curve is the exponen- 
tial damping approximation. In the right panels, black, blue, and 
red curves show the evolution of stellar obliquity, planetary, and 
stellar spins, respectively. In top panels, we assume the same ini- 
tial conditions as LWC09, with Q' g = Q', q ~ ^^^- ^^^ evolution 
for all parameters but eccentricity looks similar to LWC09's re- 
sults. The eccentricity evolution follows the exponential damping 
approximation. In middle panels, we use the same initial condi- 
tions as in the top panels, but assume a less efficient tidal damping 
in the planet {Q' g = 10^). In this case, the eccentricity evolution 
resembles the one in LWC09, but the planetary spin-orbit syn- 
chronization occurs on a longer timescale, as expected. In bottom 
panels, we use the same initial conditions as in the top panels, but 
artificially multiply the eccentricity evolution contribution from the 
planet by some small factor. Here, we recover the results of LWC09. 



artificially reducing the planetary contribution term in 
the eccentricity evolution equation (i.e., second term in 
Eq. [5]), by a factor of 10"^, i.e., completely suppressing 
the eccentricity damping effect due to the planet ^. As 
already mentioned, this discrepancy occurs because the 
numerical results in LWC09 indeed underestimated the 
eccentricity damping in the planet, since their code in- 
correctly had an additional factor of n multiplied in the 
second term of Eq. [8] (B. Levrard, private communica- 
tion). 

4.3. Lifetimes of Transiting Planets on Eccentric Orbits 

In Section 4.1, we identify two characteristic evolution- 
ary paths for Darwin-unstable planets. We now study 
the tidal evolution of eccentric transiting planets by in- 
tegrating the tidal equations forward in time with various 
tidal quality factors, and further investigate the implica- 
tions of these two paths. 

We use the currently observed parameters as initial 
conditions (see Table [1}. The orbital and stellar spin pa- 
rameters are taken from The Extrasolar Planets Ency- 
clopedia (http://exoplanet.eu/j), and references therein. 



^ Note that the choice of the factor is rather arbitrary. The 
exact number can be anything as long as the planetary contribution 
becomes negligible compared to the stellar contribution. 



For the planetary spin period we assume a small ini- 
tial value (0.2 days), but the overall results are not af- 
fected by the exact choice of this value. This is because 
the planetary spin carries a very small angular momen- 
tum, and thus pseudo-synchronization with the orbit is 
achieved very quickly (see also LWC09). For the stellar 
obliquity, we use the observed projected value A when 
RM measurements are available. Here, we implicitly as- 
sume that both the angle between the stellar spin axis 
and the line of sight (i,) and the angle between the or- 
bital angular momentum and the line of sight (io) are 
c± 90°. Thus, from cos e,. = sin?* cos A sin lo -I- cos u cosio 
(jFabrvckv fc Wiimll2009D . the stellar obliquity becomes 
comparable to the projected one (e* ~ A). For systems 
without RM measurements, we assume an initial stellar 
obliquity e*.o — 2°. This choice is rather arbitrary, but 
is motivated by the typical projected obliquity observed 
(see Table [2]). 

First, we show the tidal evolution for typical tidal 
quality factors (Q'^ g = Q'p o = 10®), with the scaling 
of Q' = QqUo/u, and without magnetic braking effect. 
Our goal here is to show how different the evolution 
can be within observational uncertainties. We exclude 
WASP-12 and HAT-P-1 from our analysis because of the 
uncertainty in their age. The evolution of semimajor 
axis and eccentricity of each system is shown in Fig. [B) 
In each panel, we show the results of 5 different inte- 
grations. Solid curves represent the results with fidu- 
cial values of a and e, while dotted curves correspond 
to 4 different combinations of maximum and minimum 
a and e values within their error bars. As expected 
from Section 4.1, the Darwin-stable systems CoRoT-3 
and HD 80606 migrate outward and inward, respectively, 
and arrive at their tidal equilibrium, with orbital decay 
eventually stopping. The borderline systems HAT-P-2 
and WASP-10 are likely Darwin-unstable, but may ar- 
rive at a stable tidal equilibrium within the observed un- 
certainties. Thus, it is important to know the orbital 
and spin parameters as well as possible to determine the 
final fate of these borderline systems. The other sys- 
tems are definitely Darwin-unstable within the current 
observational accuracy, and their planets spiral toward 
the central star on different timescales. Vertical lines 
indicate the estimated ages with uncertainties. With 
Q'* ~ QpO ~ ^^^ some systems undergo orbital de- 
cay too quickly to be compatible with their likely age 
(e.g., WASP-18, XO-3). Therefore, these results clearly 
imply that a single set of values for the tidal quality fac- 
tors cann ot reasonably apply to a ll systems (also see, for 
example, iMatsumura et al.ll2008D . 

We repeated the above calculations for various ini- 
tial stellar and planetary tidal quality factors ranging 
over the interval 10^ < Qo < 10^, with the scaling of 
Q' — Q'ono/n. In Fig. [71 we show the range of values 
for the tidal quality factors that allow planets to sur- 
vive (a > i?*), and stay eccentric (e > 0.001) for 0.1, 1, 
and 10 Gyr. Here magnetic braking is not included. For 
Darwin-stable systems (CoRoT-3, and HD 80606 in our 
list), the minimum stellar and planetary tidal quality fac- 
tors correspond to the orbital circularization times. For 
all Darwin-unstable transiting systems with noncircular 
orbits, we find the same trend as in Section 4.1: the circu- 
larization time is largely determined by the dissipation 



in the planet, while the survival time is largely deter- 
mined by dissipation inside the star. In other words, for 
Darwin-unstable planets, the minimum planetary tidal 
quality factors can be inferred from the circularization 
time, while the minimum stellar tidal quality factors can 
be inferred from the orbital decay time. We demonstrate 
this below. 

The approximate minimum planetary tidal quality fac- 
tor that allows a planet to keep a non-circular orbit (e > 
0.001) for a certain time (in our examples, 0.1 — 10 Gyr) 
can be determined by assuming that the eccentricity evo- 
lution depends only on tidal dissipation inside the planet 
(de/dt) ~ {de/dt)p. We rewrite Eq. [S] as follows by 
assuming pseudo-synchronization of the planetary spin 
{dujp/dt = 0), as well as conservation of angular momen- 
tum (a(l — e^) = const): 
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By integrating the above equation from the currently ob- 
served eccentricity to e = 0.001, and solving for Q' q, we 
obtain the minimum planetary tidal quality factors to 
keep a non-circular orbit for 0.1, 1, and 10 Gyr, respec- 
tively. These values are plotted as the vertical dashed 
lines in Fig. [71 

Similarly, we can determine the approximate minimum 
stellar tidal quality factor that allows a planet to survive 
for a certain time before plunging into the central star, by 
assuming that the semimajor axis evolution depends only 
on tidal dissipation inside the star (da/dt) ~ {da/dt)„. 
Note that this is a reasonable approximation when the 
pseudo-synchronization of the planetary spin is achieved, 
and the orbit is nearly circular. We rewrite Eq. [71 as 
follows by setting e = 0, 
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We integrate the above equation from a = ao(l — Bq) to 
i?* , and solve for Q'^ q to obtain the horizontal dashed 
lines. Here, we make use of the fact that the difference 
in eccentricity damping times does not strongly affect 
the orbital decay time, and assume that the orbital de- 
cay time of any eccentric Darwin-unstable system can be 
well described by that of a system with equal angular 
momentum and a circular orbit. 

As seen in Fig. [71 the agreement of both horizontal and 
vertical dashed lines with the calculations for Darwin- 
unstable systems is very good. Since Darwin-stable sys- 
tems (CoRoT-3 and HD 80606) arrive at the stable tidal 
equilibria and stop migrating, the horizontal dashed lines 
for these systems significantly differ from the calculated 
results. Now we present some example runs along hor- 
izontal and vertical lines for WASP- 17 to better under- 
stand their implications. The left panel of Fig. [HI shows 
three runs along the lowermost horizontal line that corre- 
sponds to the survival time of 0.1 Gyr. Dotted, solid, and 
dashed curves show the evolutions with the same initial 
stellar tidal quality factor Q'^ q = 9.13 x 10^, and different 
initial planetary tidal quality factors QLq — 7.43 x 10^, 
7.43 X 10^, and 7.43 x 10'^, respectively. The figure con- 
firms that the orbital decay time is largely determined 
by the tidal dissipation in the star, since there is no 
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obvious difference depending on Q' g values. At the 
vertex of vertical and horizontal dashed lines in Fig. [7] 
(i.e., (gt^o, Q;^o) = (9.13 X 10*, 7.43 x 10^)), we find 
that the semimajor axis and eccentricity damp roughly 
on the similar timescale {t^ ^ Ta ^ 0.1 Gyr). For the 
smaller Q' q, the eccentricity damps much faster than 
the orbital decay {t^ < Ta ^ 0.1 Gyr), while for the larger 
Q'pQ, the eccentricity damps slower than it is expected 
from the other two curves, and on a similar timescale to 
the orbital decay (te ^ Tq ^ 0.1 Gyr). Similarly, the right 
panel of Fig. |8] shows three runs with the same initial 
planetary tidal quality factor Q'^q = 7.43 x 10^, and dif- 
ferent initial stellar tidal quality factors Q'^ q = 9.13x 10^, 
9.13 X 10"*, and 9.13 x lOl When Q'^^ is smaller than the 
vertex value, the orbit decays much faster than 0.1 Gyr, 
and the circularization happens on the similar timescale 
(te '^ Ta < 0.1 Gyr). On the other hand, when Q'^ q is 
larger than the vertex value, the orbit decays much slower 
than 0.1 Gyr, and the circularization time is shorter than 
the orbital decay time (rg ~ 0.1Gyr< r^). The fig- 
ure also implies that the orbital circularization time is 
largely determined by the dissipation in the planet unless 
Te ^ Ta- Thus, wc find that Tg ~ Tq is a good approxima- 
tion along the horizontal dashed lines, while r^ < r^ is 
a good approximation along the vertical lines. In other 
words, the region below the diagonal line drawn by con- 
necting the vertices of vertical and horizontal lines is 
the stellar-dissipation-dominated region, while the region 
above the line is the planetary-dissipation-dominated re- 
gion. 

We repeat these integrations by including magnetic 
braking effects. As can be seen in Fig. [9l there is lit- 
tle difference between the cases with and without mag- 
netic braking. Note that dashed lines are the same as the 
ones in Fig. [71 and th us do not take account o f the mag- 
netic braking effects. [Barker fc Ogilvia l|2009l ) suggested 
that the effect of magnetic braking can be important for 
systems with rapidly spinning stars (w+cose/n 3> 1). 
From Table [21 there are five such systems (CoRoT-3, 
CoRoT-6, HAT-P-2, HD 80606, and WASP-7). Among 
them, CoRoT-3, CoRoT-6, HD 80606, and WASP-7 are 
Darwin-stable systems, while HAT-P-2 is a borderline 
case that can be either Darwin-stable or unstable within 
observational uncertainties. Out of these systems with 
fast spinning stars, CoRoT-3, HAT-P-2, and HD 80606 
have eccentric planets and are shown in Fig.[9l Excluding 
Darwin-stable cases (CoRoT-3 and HD 80606), HAT-P-2 
indeed shows a significant difference between Fig. [71 and 
[9lfor the 1 and 10 Gyr cases. 

5. PAST TIDAL EVOLUTION OF TRANSITING PLANETS 

5.1. Observational Implications on the Origins of 
Close-in Planets 

In Section 1, we mentioned the two main scenarios to 
form close-in planets — planet migration in a disk, and 
tidal circularization of an eccentric orbit, which may be 
obtained as a result of planet-planet scattering or Kozai- 
type perturbations. It is nontrivial to determine which 
formation mechanism dominates, but there are at least 
a few observational indications that support the second 
scenario. 

One of them relates to the orbita l distribution of plan- 
etary systems. [Wright et al.l ([2009[ ) compared the prop- 
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Fig. 6. — The evolution of semimajor axis (black curves) and 
eccentricity (orange curves) for Q'^ q = Q\ 



p,0 



lO*". Tidal qual- 
ity factors scale as Q' = QgUo/n, and magnetic braking is not 
included. Solid curves correspond to the fiducial values, while dot- 
ted curves correspond to four different combinations of maximum 
and minimum semimajor axis and eccentricity, allowed within the 
uncertainties. Vertical lines show the age of each system (dashed 
lines) with uncertainties (dotted lines). As expected from Sec- 
tion 4.1, CoRoT-3 and HD 80606 arrive at their stable tidal equi- 
librium, while the other planets spiral into the central star. This 
figure clearly demonstrates that different tidal quality factors ap- 
ply to different systems, since some planets fall within the Roche 
limit of their stars on time scales much shorter than the age of the 
star with a common value of Q'^ g. 
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Fig. 7. — Combinations of stellar and planetary tidal quality 
factors which keep a non-zero eccentricity and allow survival of 
the planet in forward integration of the tidal equations for 0.1, 1, 
and 10 Gyr (orange, blue, and red regions, respectively). Magnetic 
braking is not included, and tidal quality factors change as Q' = 
QqUo/u. Vertical, and horizontal dashed lines are determined by 
assuming (de/dt) ~ {de/dt)p, and (da/dt) ^ {da/dt)^, respectively 
(see text for details). The system's lifetime is largely determined 
by the tidal dissipation in the star, and the circularization by that 
in the planet. 



erties of multiple-planet systems with single-planet sys- 
tems (i.e., with no obvious additional giant planet), and 
showed that their semi-major axis distributions are dif- 
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Fig. 8. — Tidal evolution of WASP-17 with different initial tidal 
quality factors along the dashed lines in Fig. [7] Black and orange 
curves correspond to semimajor axis and eccentricity evolutions, 
respectively. The vertical dashed lines are drawn at 0.1 Gyr for 
comparison. Left: Different initial conditions along the lowermost 
horizontal dashed line in Fig. [7] that indicates the survival time 
of 0.1 Gyr. Three runs with the same initial stellar tidal quality 
factor Q'^ q = 9.13x10'*, and different initial planetary tidal quality 
factors are shown. Dotted, solid, and dashed curves correspond 
to Q'pQ = 7.43 X 10*, 7.43 x 10^^, and 7.43 x 10^, respectively. 
Orbital decay time is determined by the tidal dissipation in the 
star, since the decay time does not change for different Q' q values. 

For Q' Q > 7.43 x 10^, it is clear that the eccentricity damps on 
a similar timescale to the orbital decay. Right: Different initial 
conditions along the leftmost vertical dashed line in Fig. [7] that 
indicates the circularization time of 0.1 Gyr. Three runs with the 
same initial planetary tidal quality factor Q' q = 7.43 X 10^, and 
different initial stellar tidal quality factors are shown. Dotted, 
solid, and dashed curves correspond to Q^ q = 9.13x10'^, 9.13x10*, 

and 9.13x10^, respectively. For Q'^q < 9.13x 10*, both the orbital 
decay and the circularization times are much shorter than 0.1 Gyr. 
For Q'^ g > 9.13 x 10*, both become comparable to, or longer than 
0.1 Gyr. 
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Fig. 9. — Same as Fig.fT] but with the effects of magnetic braking 
included. There is very little difference for the future tidal evolu- 
tions with and without magnetic braking. Vertical and horizontal 
dashed lines are the same as in Fig. [7] 



ferent. While single-planet systems have a double-peaked 
distribution, which is characterized by a pile-up of giant 
planets between 0.03 and 0.07 AU (the so-called "3-day 
pileup") as well as a jump in the number of planets be- 
yond 1 AU, multiple-planet systems have a much more 
uniform distribution. They also pointed out that the 
occurrence of close-in (a < 0.07 AU) planets is lower for 



multiple-planet systems, and that planets beyond 0.1 AU 
in multiple-planet systems exhibit somewhat smaller ec- 
centricities compared to the corresponding single ones. If 
confirmed by future observations, this trend would favor 
planet-planet interaction scenarios over a migration one, 
because there is no obvious reason why the orbital dis- 
tributions of single- and multiple-planet systems should 
be different for planet migration. On the other hand, 
gravitational interactions combined with tidal circular- 
ization may be able to explain such a difference, be- 
cause strong gravitational interactions tend to disrupt 
the system, and thus currently observed multiple-planet 
systems are unlikely to have been strongly perturbed by 
stellar/planetary companions, and/or to have undergone 
catastrophic scattering events. 

Another indication comes from the similarity between 
the stellar obliquity distribution derived from the ob- 
served systems, and the distributio n predicted by the 
Kozai migratio n scenario JF abrvck v fc Tremaind 120071 ; 
IWu et al.ll200l . iTriaud et all (|2010i) observed the RM 
effect for six transiting hot Jupiters, and found that four 
of their targets appear to be in retrograde orbits with 
a projected stellar obliquity > 90°. Combining the pre- 
vious 20 systems with such measurements, they pointed 
out that 8 out of 26 systems are clearly misaligned, and 
that 5 out of 8 misaligned systems exhibit retrograde 
orbits. They also derived the stellar obliquity distribu- 
tion by assuming an isotropic distribution of the stellar 
spin with respect to the line of sight, and found that 
the distribution closely ma tches that e xpected from the 
Kozai mi g ration sc enario (iFabrvcky fc Tre mainc 2003; 



Wu et al. 



Wu et al. 



200J. IFabrvcky fc Tremaind (|2007l ) and 
(I2007f ) independently studied the possibility of 



forming a close-in planet by considering the combined 
effects of secular perturbations due to a highly-inclined 
distant companion star (i.e., Kozai-type perturbations) 
and tidal interactions with the central star. In this sce- 
nario, a Jupiter-mass planet which was initially on a 
nearly-circular orbit at ^ 5 AU can become a hot Jupiter. 
The mechanism involves a companion star on a highly- 
inclined orbit (~ 90°), which perturbs the planetary or- 
bit and increases its eccentricity. Once the pericenter 
distance of the planet becomes small enough for tidal in- 
teractions with the central star to be important, energy 
dissipation leads to circularization of the planetary or- 
bit, and eventually to formation of a hot Jupiter with a 
small, or nearly zero eccentricity. They found that hot 
Jupiters formed this way tend to be in misaligned or- 
bits, and frequently in retrograde orbits. Planet-planet 
interactions around a single star (without a binary com- 
panion) could also form hot Jupiters via Kozai migration 
(jNagasawa et al.l[200 8). 

Yet another clue regarding the origin of close-in plan- 
ets is related to the above scenario, and comes from the 
inner edge of the orbital distribution for hot Jupiters. 
iFord fc Rasiol ([2006) proposed that the observed inner 
cutoff for hot Jupiters is defined not by an orbital pe- 
riod, but by a tidal limit, and studied such a cutoff of the 
distribution of close-in giant planets in the mass-period 
diagram by performing a Bayesian analysis. Assuming 
that a slope in such a diagram follows the Roche limit, 
they found that the observations suggest an inner cutoff 
close to twice the Roche limit. This can be explained 
naturally if the planetary orbits were initially highly ec- 
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Fig. 10. — Planetary and stellar mass ratio as a function of 
semi-major axis normalized to the Roche-limit distance. For non- 
transiting planets (orange circles) without planetary radius infor- 
mation, we assume a Jupiter radius for planets with mass larger 
than 0.1 Mj, and a Neptune radius for planets with smaller mass. 
Most close-in planets lie beyond 2 times their Roche Limit. One 
planet, WASP-12b, has a non-zero eccentricity, and a very small 
semi-major axis (a ~ l.Sajj). 



centric, and later circularized via tidal dissipation while 
conserving orbital angular momentum. They suggested 
that this result is inconsistent with the migration sce- 
nario, because migration would lead to an inner edge 
right at the Roche limit (a factor of 2 further in than 
what is observed). 

In Fig. [ini we extend the work of iFord fc Rasiol l|2006[) 
by including all more recently discovered planets, and we 
plot planetary and stellar mass ratio against s emi-major 
axis in terms of the Roche limit separation (jPaczvhskil 

[Unl) 



Or 



:Rv 



M„ 



3 (M* + Mp) 



-1/3 



Rp fMp 



0.462 \M,. 



-1/3 



(19) 

Here, we assume that the Roche radius of the planet, 
which is defined so that its spherical volume is equal 
to the volume within the Roche lobe, is equal to the 
planetary radius. Thus, the Roche limit separation used 
here is the lower limit. For non-transiting planets with- 
out a measured planetary radius (orange circles), we as- 
sume a Jupiter radius for planets with mass larger than 
0.1 Mj, and a Neptune radius for less massive planets. 
It is obvious that most planets still exist beyond twice 
the Roche limit (vertical dashed line). However, there 
are now 5 planets which lie within this limit (OGLE- 
TR-56, CoRoT-1, WASP-4, WASP-19, and WASP-12) 
with a/a-R, ^ 1.70, 1.67, 1.66, 1.30, and 1.24, respectively. 
We need to assess whe t her th ese planets will change the 
claim bv lFord &: Rasiol (|2006l ). and support the migration 
scenario over the scattering/Kozai-cycle scenario. Note 
that the two recently discovered "extreme" close-in plan- 
ets, with orbital periods less than 1 day, CoRoT-7 b and 
WASP-18 b, have a/aji ^ 2.76, and 3.52, respectively. 
The existence of at least some of these planets inside 



2afl may still be explained as a result of tidal circular- 
ization of an eccentric orbit. One of the possibilities is 
that the orbits of these planets were originally circular- 
ized beyond twice the Roche limit, but the planets have 
migrated inward after circularization due to tidal dissi- 
pation in the star. All of the planets within 2a/j can 
be potentially explained this way. Another possibility 
is that their orbits used to have a pericenter distance 
close to the Roche limit, but the initial eccentricity was 
smaller, e ~ 0.7. In such a case, the expected period of 
the circularized orbit would be comparable to the current 
orbital period for OGLE-TR-56, CoRoT-1, and WASP- 
4, since they all have a/oR ~ 1.7. However, systems with 
smaller ratios of semimajor axis to Roche limit separation 
(WASP-12 and WASP-19) are unlikely to be explained 
this way. This is because their current semimajor axes 
(a/oR ~ 1.24 and 1.30) would demand small initial ec- 
centricities (ci ^ 0.24 and 0.30), and thus small initial 
semi-major axes {oq ~ 0.024 AU and 0.018 AU, respec- 
tively). This means that they have to be either born at 
such close-in locations, or scattered into such an orbit, 
which would be very difficult. Yet another possibility is 
that these planets used to have a smaller Roche limit 
separation, due to a larger mass, or a smaller radius. 
The orbits of these planets might have been circular- 
ized as the planet lost ma ss (e.g. iLammer et al. 200^ 
iLecavelier Pes Etang£ll2007l ). or the planetary radius ex- 
pand e d due to tidal heating (e.g. iBodenheimer et al] 
[200H iGu et all [2003h . For OGLE-TR-56, CoRoT-1, 
WASP-4, WASP-19, and WASP-12 to be circularized at 
twice the Roche limit, either the past planetary masses 
must have been 2.13, 1.78, 1.53, 4.19, and 5.98 Mj, re- 
spectively, or the past planetary radii must have been 
1.18, 1.24, 1.21, 0.829, and 1.11 i?j, respectively. Since 
the mass- loss rate can only be at most ^ 10%, even for a 
low-density planet like WASP-12 (jLammer et al. 2009), 
it is unlikely that a larger mass in the past could be the 
correct explanation. On the other hand, tidal infiation 
of the planetary radius is a transient phenomenon (e.g. 
Ilbgui fc Burrowa l2009[ ) . To explain the current orbital 
radius by invoking a smaller planetary radius in the past, 
we have to catch the planet just as its radius is inflat- 
ing. Such a scenario may be possible, but appears un- 
likely. Interestingly, CoRoT-1 is obser ved to be strongl y 
misaligned with A = — 77± 11 degrees (jPont et al.ll2010f ). 
which suggests a violent origi n, while WASP-4 ha s a stel- 
lar obliquity of 4^43 degrees (jTriaud et al.l 120101 ) . which 
is consistent with alignment. We urge observers to de- 
termine the projected stellar obliquity for OGLE-TR-56, 
WASP-19, and WASP-12. Alignment does not necessar- 
ily support the planet migration scenario over the violent 
origin, but misalignment would clearly imply a significant 
past dynamical interaction involving scattering or Kozai 
cycles. 

5.2. Past Evolution of Transiting Planets 

As we pointed out in the previous subsection, planet- 
planet (or planet-companion star) interactions may well 
be responsible for the majority of close-in planets. Now 
we further explore this possibility by performing inte- 
grations of the tidal evolution equations backward in 
time. Particular focus is on the differences in evolution 
of Darwin-unstable planets depending on two evolution 
paths. 
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iJackson et al.1 (|2008f ) first performed such a study, and 
estimated the "initial" eccentricity distribution of close- 
in planets. They found that this distribution agrees well 
with the observed eccentricity distribution of exoplanets 
on wider orbits, and they proposed that most currently 
observed close-in planets had considerably wider, and 
more eccentric orbits in the past. However, their study 
had many uncertainties. First, the possible effects of 
other planets and binary companions were neglected. As 
we pointed out in Section 2, these effects are most likely 
unimportant at present, but they could well have been 
playing a crucial role in the past, especially if another 
object was responsible for Kozai-type perturbations, or 
gravitational scattering. Second, the evolution of stars 
and planets is neglected. For example, both planetary 
and stellar radii may have been significantly different in 
the past. Moreover, many stars could have lost a large 
amount of their spin angular momentum via magnetic 
braking. Finally, such backward dynamical integrations 
of evolution equations with energy dissipation are known 
to be diverging, and thus small changes in initial condi- 
tions can lead to much larger changes in the calculated 
"initial" values. 

With these caveats in mind, we now re-examine the 
possible past histories of transiting planets. Since Eq. [T]- 
111! are time-invariant, we can study the past evolu- 
tion by integrating the differential equations "backward 
in time" (i.e., by taking the negative of these differ- 
ential equations and integrating them from t ~ to 
Tagc). For simplicity, we assume that the planetary 
spin is pseudo-synchronized with the orbit at all times 
(i.e., dujp/dt = 0). First, we show typical results for 
WASP-18 in Fig.[TTJ Here, the initial tidal quality factors 
are Q^ = QpQ = 10^ with the scaling of Q' = QqUo/ti, 
and the magnetic braking effect is neglected. As in Fig. [51 
solid curves show the results of backward evolution with 
the fiducial orbital parameters, while four dotted curves 
correspond to different combinations of maximum and 
minimum semimajor axis and eccentricity within uncer- 
tainties. Vertical dashed and dotted lines correspond 
to the estimated stellar age with uncertainties. As ex- 
pected, both semimajor axis and eccentricity values differ 
significantly within uncertainties at the zero age (vertical 
dashed line) , or at any specific time within the age uncer- 
tainties (between vertical dotted lines). However, we find 
that the fiducial case still provides a representative evolu- 
tionary path within the stellar age uncertainties. Fig. [12] 
presents similar backward evolutions for all the systems 
shown in Fig. [5) It is clear that, except the Darwin-stable 
CoRoT-3, WASP-18 has the largest spread in orbital pa- 
rameters within uncertainties. The backward evolution 
of the other systems is largely independent of the exact 
values of orbital parameters. 

Now we repeat backward integrations of the tidal equa- 
tions by adopting various combinations of initial stellar 
and planetary tidal quality factors ranging over 10^ < 
Q'o < 10^. This allows us to study the "initial' orbital 
properties at the zero stellar age {rage = 0) within the 
estimated age uncertainties. As in Section 4.3, we also 
assumed that each system initially has the currently ob- 
served orbital and rotational parameters. 

We select systems that have solutions such that e < 1 
somewhere in the interval Tage,rmn <t< Tage.max, and 
we plot the corresponding Q'^ q and Q' q values in Fig. 1131 



Again, the tidal quality factors change as Q' — QqUq/tl^ 
and magnetic braking is neglected. Blue, green, or- 
ange, and red regions represent the maximum reachable 
"zero- age" semi- major axis for each system being 0.1, 1, 
10 AU, and > 10 AU, respectively. Vertical and horizon- 
tal dashed lines are the same as in Fig. [71 and indicate 
the minimum Q' q and Q'^ q required to have circular- 
ization and orbital decay times (i.e., the future survival 
time) of 0.1, 1, and lOGyr, respectively. As expected, 
when tidal dissipation is inefficient in both planet and 
star (top right corner area of the figure) , the planets are 
generally expected to have stayed where they are now for 
a long period of time. Lack of a significant change in its 
orbit, especially eccentricity, may be consistent with the 
initial orbital properties expected from planet migration 
scenario. For the future survival time much longer than 
10 Gyr, only such a solution with little migration is al- 
lowed in the parameter space. On the other hand, when 
tidal dissipation is highly efficient (bottom left corner), 
a planet could have started on a wide, eccentric orbit, 
which was then circularized over the lifetime of the sys- 
tem. The property of this area is in good agreement with 
the expectation from Kozai cycles and/or planet-planet 
interactions followed by tidal dissipation. The compar- 
ison of these regions with dashed lines imply that most 
planetary systems have a wide range of allowed parame- 
ter space for (Q^ g, Q' g) which is consistent with having 
started on a wide, eccentric orbit, and with comfortable 
future survival times ^1 — 10 Gyr. Th us, our results 
agree with the suggestion first made by IJackson et al.l 
(2008) ■ and show that there is a broad parameter space 
which supports the tidal circularization scenario as the 
dominant origin of close-in planets. 

When we demand that the future survival time must 
be comparable to, or slightly longer than the estimated 
stellar age. Fig. [T31 implies an interesting trend. In the 
stellar-dissipation dominated region, planets tend to have 
an initial orbit which is similar to the current one, while 
in the planetary-dissipation dominated region, planets 
tend to have an initially highly eccentric and wide orbit. 
This is clear by comparing these two regions (below and 
above the diagonal line that can be drawn by connect- 
ing the vertices of vertical and horizontal dashed lines) 
at around their survival times. Most systems in the fig- 
ure except WASP-10, WASP-14, WASP-18, and WASP-6 
have the estimated age of ^ 1 — 10 Gyr, while WASP- 
10, WASP-14, and WASP-18 have ~ 0.1 - 1 Gyr, and 
WASP-6 has > 10 Gyr. Investigating the corresponding 
regions, we find that the maximum zero-age semi-major 
axis can be very large (upto > 10 AU) in the planetary- 
dissipation dominated region. On the other hand, such a 
solution is less likely in the stellar-dissipation dominated 
region, and the area with little change in orbital radii 
tends to occupy the largest parameter space. 

Fig. [14] presents similar results with magnetic braking 
effects included. As we can see, some systems are much 
less affected by magnetic braking than others. Our re- 
sults show that, when either the stellar spin is slow, or 
the mass ratio is very low, the evolution is less affected 
by magnetic braking. Clearly affected systems include 
CoRoT-3, HAT-P-2, HD 189733, WASP-10, WASP-14, 
WASP-18, WASP-5, and XO-3, for which the average 
stellar spin period is ~ 8 days and the average mass ra- 
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Fig. 11. — The backward evolution of semi-major axis a 
(black curves) and eccentricity e (orange curves) for WASP-18 with 
Q't = Q'p ~ W^. Solid curves are the nominal values, and dot- 
ted curves show four independent runs with different combinations 
of a and e within uncertainties. Vertical lines show the estimated 
stellar age with uncertainties. 
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Fig. 13. — Combinations of stellar and planetary tidal quality 
factors that allow a planet to survive (i.e., e < 1) for the stellar age 
with uncertainties in backward integration of the tidal equations. 
Magnetic braking is not included, and tidal quality factors change 
as Q' = Qgno/n. Blue, green, orange, and red areas correspond to 
a maximum zero-age semi-major axis of 0.1, 1, 10, and > 10 AU, 
respectively. Also plotted are the same vertical and horizontal 
dashed lines shown in Fig. [7] 
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Fig. 12. — Same as the previous figure for other transiting sys- 
tems. The evolution of semi-major axis and eccentricity is shown 

with Q',0 = q;,o = lo''. 



tio is ~ 0.006, while the corresponding values are about 
36 days and 0.001 for the others. 

5.3. Stellar Obliquity Evolution 

Fig-HH and [TBI show similar plots to Fig.[T31 and [HI for 
the stellar obliquity. Blue, green, orange, and red areas 
correspond to a maximum possible zero-age obliquity of 
5°, 20°, 40°, or > 40°. Again, also plotted are the same 
vertical and horizontal lines as in Fig. jT] 

In Fig. dSl and (H CoRoT-3, HAT-P-2, HD 17156, 
HD 80606, WASP-14, WASP-17, WASP-18, WASP- 
5, WASP-6, and XO-3 have RM measurements, and 
thus known projected stellar obliquities. Among them, 
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Fig. 14. — Same as Fig. 1131 but with the effects of magnetic 
braking included. Systems that are affected the most by magnetic 
braking have the rapidly rotating star, and/or a relatively high 
mass ratio. For CoRoT-3, we don't get any solutions. 



CoRoT-3, HD 80606, WASP-14, WASP-17, and XO-3 

have a significant misalignment (> 20°). Naturally, the 
maximum possible misalignment at zero age can be very 
large (> 40°) for these systems with almost any values 
of (Q'*i Q'p)- Oil the other hand, for systems with small 
measured projected obliquities (HAT-P-2, HD 17156, 
WASP-18, WASP-5, and WASP-6), or systems with no 
RM measurements (CoRoT-5, GJ 436, HAT-P-11, HAT- 
P-13, HD 189733, and WASP-10), we find either little or 
large change in obliquity over the stellar age. 
Consider a few specific cases. For HAT-P-2 and HAT- 
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P-11, there are no solutions for obliquity larger than 5°. 
This indicates that, for HAT-P-2, the "zero-age" obliquity 
was likely similar to the current nominal value A = 1.2°. 
However, since HAT-P-2 's measured obliquity has a very 
large uncertainty (A — 1.2° ± 13.4°), it is possible that 
the actual present value is much larger. Note that even 
if the stellar obliquity turns out to be small, HAT-P-2 is 
still likely to have the scattering/Kozai-cycle origin. This 
is partly because of its high eccentricity (e ^ 0.52), and 
partly because of a large range of tidal quality factors 
that allow a much wider orbit in the past (see Fig. [T51 
and [HI). For HAT-P-11, there are no RM measurements 
so far, but the system shows a similar result to HAT-P- 
2 with e,^o ~ 2°. If the current obliquity turns out to 
be < 2°, our plot indicates that HAT-P-11 would not 
have had a large obliquity in the past. This in turn indi- 
cates that, if HAT-P-11 initially had a large stellar obliq- 
uity, we should be able to see a clearly misaligned orbit 
through future observations, independent of the actual 
tidal quality factors for the system. The orbital eccen- 
tricity of the planet is e = 0.198 ± 0.046, which could 
have been produced either via planet-disk, or planet- 
companion interactions. Interestingly, Fig. [13] and [14] 
show that the orbital radius of HAT-P-11 b is unlikely 
to have been changed significantly via tidal evolution, 
unless the tidal dissipation inside the planet is rather ef- 
ficient (Q' < 10^). If the tidal dissipation is inefficient, 
our result suggests that HAT-P-11 b is likely to have the 
migration origin. We have to wait for future observations 
to estimate whether the planet is likely formed via disk 
migration, or tidal circularization of a highly eccentric 
orbit^. 

For all the other systems, the zero- age obliquity can 
be as high as > 40° if stellar tidal quality factors are 
relatively small. More specifically, such solutions are al- 
lowed in the stellar-dissipation dominated region with 
Te ~ Ta- The stcUar obliquity generally damps on a 
similar timescale to th e orbital decay (see, e.g., LWC09, 
[Barker fc Ogilviell2009l as well as Section 4.1). However, 
this plot demonstrates that the stellar obliquity could be 
damped from high (> 40°) to low {^ 2°) values within 
the current stellar age. 

Note that, these small stellar tidal quality factors 
which allow the fast damping of stellar obliquities, also 
lead to relatively short survival times for planets. By 
comparing these zero-age high obliquity regions with the 
horizontal dashed lines, we find that these tidal qual- 
ity factors lead to survival times at most comparable to 
the current stellar age. Or in other words, relatively 
small parameter spaces are allowed for these high zero- 
age obliquities. For example, when we demand that the 
expected survival time in forward tidal evolution must 
be comparable to or larger than the stellar age, we find 
that the stellar tidal quality factor for CoRoT-5 must 
be Q'^ > 10^. This includes a small region of the red, 
large "zero-age" stellar obliquity area in Fig. [TS] How- 
ever, if Q'^ > 10^, such an area disappears. In the lat- 

® As we were writing up this draft, two groups released the resuhs 
that HAT-P-11 is a highly inclined system. The sky-projected 
stellar obliquity is A = 103_]^g degrees by I Winn et al.l l|2010d ). and 
103_jg degrees bv lHirano et al.l ||201CII '). Thus, the system is likely 
to have the scattering/Kozai-cycle origin, rather than the migration 
one. 



ter case, if the observations find a small stellar obliquity 
for CoRoT-5, it is unlikely that the stellar obliquity was 
much larger in the past. In some cases, such a region 
with a high zero-age obliquity and a comfortable future 
survival time may not even exist. For HAT-P-13, a simi- 
lar comparison shows that the stellar tidal quality factor 
must be Q'^ > 10^ for the survival time to be comparable 
to or larger than the stellar age. This corresponds to the 
blue, small "zero-age" stellar obliquity area in Fig. \TE\ 
and thus we expect that the stellar obliquity of HAT- 
P-13 could not have been changed much due to tidal 
evolution (i.e., similar to HAT-P-2 or HAT-P-11). In- 
terestingly, HAT-P-13 is the only system in our sample 
that has an additional planetary companion. Thus, the 
system may be an example of a close-in planet formed 
via migration scenario. If this is the case, we predict the 
future observations will find a small stellar obliquity for 
HAT-P-13 ^. 

In short, by comparing Fig. [13] with Fig. 1151 we find 
that the evolution history is largely divided into two 
cases, depending on the relative efficiency of energy dis- 
sipation inside the star and the planet. In the planetary- 
dissipation dominated region, the planets could have had 
a wide, eccentric orbit in the past, and the stellar obliq- 
uity damps on a similar timescale to the orbital decay. 
The initial conditions implied in this region are consistent 
with those expected from the scattering/Kozai-cycle ori- 
gin of the close-in planets. On the other hand, in the 
stellar-dissipation dominated region, the system could 
have had a large stellar obliquity in the past, although 
the initial orbital radius and eccentricity are likely sim- 
ilar to the current values, as planet migration scenario 
would suggest. The further inspection of these figures 
show that a unique evolution is possible in the transition 
region of these two, where the dissipation effects due to 
a planet and a star are comparable. There, the system 
could have had a wide, eccentric orbit, as well as a large 
stellar obliquity in the past. Our results also imply that 
if the currently observed stellar obliquity distribution 
is due to Kozai migration as suggested by iTriaud et al.l 
(2010), then most exoplanetary systems have planetary- 
dissipation dominated tidal interactions (i.e., Tg < t^). If 
Te ~ Ta for most planetary systems, the current obliquity 
distribution should be very different from the initial one, 
and the dynamical history prior to the tidal dissipation 
should be wiped out. 

Recently, IWinn et al.l (|2010a[ ) suggested that the stel- 
lar obliquity is preferentially large for hot stars with ef- 
fective temperatures Tcs > 6250 K, and proposed that 
such a trend can be explained because photospheres of 
cool stars can realign with the orbits due to tidal dis- 
sipation in the convective zones, without affecting the 
orbital decay. Our results suggest that yet another pos- 
sibility might be that the evolution of the systems with 
hot and cool stars correspond to slow and fast obliquity 
damping regions, respectively. If that is the case, the 
stellar obliquity may stay similar to the original value 
for the planetary system with a hot star, because tidal 
dissipation is dominated by the planet (i.e., Te < r^), 
while the obliquity may be damped to a small value for 
the system with a cool star, because tidal dissipation 

^ Indeed, a recent observation found a well-aligned orb it for 
HAT-P-13b with A = -1.9 ± 8.6 degrees l lWinn et al.ll2010bl) . 
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Fig. 15. — Combinations of stellar and planetary tidal quality 
factors that allow a planet to survive (i.e., e < 1) for the stellar 
age with uncertainties in backward integrations of the tidal equa- 
tions. Magnetic braking is not included, and tidal quality factors 
are scaled as Q' = Q'gno/n. Blue, green, orange, and red areas 
correspond to a maximum zero-age stellar obliquity of 5, 20, 40, 
and > 40 degrees, respectively. Also plotted for reference are the 
same vertical and horizontal dashed lines shown in Fig. [7] 
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Fig. 16. — Same as Fig. 1151 but with magnetic braking included. 
Systems which are affected the most have the rapidly rotating star, 
and/or a relatively high mass ratio. For CoRoT-3, we don't get any 

solutions. 



is either dominated by the star (i.e., Tg ~ Tq), or the 
star and the planet have the comparable dissipation ef- 
fects. In our samples of eccentric systems, CoRoT-5, 
GJ 436, HAT-P-11, HAT-P-13, HD 17156, HD 189733, 
HD 80606, WASP-10, WASP-5, and WASP-6 have the 
effective temperatures less than 6250 K. For the range of 
tidal quality factors we use, we don't see any trend that 
these cool systems prefer t^ ^ Ta- One possibly interest- 
ing example is HD 17156, for which the large parameter 
space is allowed for the stellar-dissipation dominated case 
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Fig. 17. — Same as Fig. [7] but with different scalings for tidal 
quality factors. 



6. DIFFERENT SCALES FOR TIDAL QUALITY FACTORS 

So far, we have focused on a scaling of Q ex 1/n. 
However, this scaling is not appropriate unless a close- 
in planet system reaches a true dual synchronization of 
n ~ u^ ^ ujp. This cannot be reached unless both ec- 
centricity and obliquities are zero, or some extra torques 
are acting on the system. For a planetary spin, although 
the true synchronization with the orbit normal does not 
happen until the final spiral-in of the planet, the pseudo- 
synchronization reaches quickly (see e.g., Section 4.1 and 
4.2). On the other hand, a stellar spin changes on a sim- 
ilar timescale to the orbital decay, and thus it's reason- 
able, for most systems, to assume that the stellar spin 
is far from synchronization. As we have seen in Sec- 
tion 3.1, the semi-diurnal tide with the forcing frequency 
of 120; — 2n I dominates the energy dissipation before the 
spin-orbit synchronization, while the annual tide with 
|2w — n\ takes over once the synchronization approaches. 
Therefore, in this section, we scale the planetary and 
stellar tidal quality factors as Qp oc l/\2ujp — n\ and 
Q* oc l/|2w* — 2n|, respectively, and investigate the dif- 
ferences in future and past evolution from the results of 
the scaling Q ex 1/n. 

In Fig. 1171 we repeat the similar tidal evolution forward 
in time as Fig. [7| with the initial tidal quality factors 



{Te ^ Ta). 



ranging over the interval 10'' < Q'q < 10^. The figure has 
a similar general trend to Fig. [71 but the larger parameter 
space is allowed. The vertical and horizontal lines are 
obtained by rescaling Eq. 1171 and fTSl bv no/\2ujpfl — no\ 
and no/|2w*^o ~ "-ol: respectively. The agreement with 
the integration of the complete tidal equations is pretty 
good. 

Similarly, Fig. (TB] and \W\ show the corresponding re- 
sults to Fig. [13] and [151 respectively, with the modified 
Q scalings. Both of these figures have a good agreement 
with the Q oc 1/n case. 

7. DISCUSSION AND CONCLUSIONS 

Close-in planets may be formed via planet migration 
in a disk, or via tidal circularization of a highly eccentric 
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Fig. 18. — Same as Fig. 1131 but with different scalings for tidal 

quality factors. 



,-. CO 

60 CD 
O 

^ LD 



00 CD 
O 



^^^^ 


WKm.-^%^--. 


RH 


1^- p - 1 3i|||||||||H|H^^ 




|^^^^^Ej^^^^EF|T7771 






||||jj^-- = X0-3|||| 



bC CD 

— lO 



456789i56789t567B»156789 
log(Q'^) log(Q;) log(Qg log(Q'p) 

Fig. 19. — Same as Fig. 1151 but with different scalings for tidal 
quality factors. 



orbit following planet-planet scattering or other gravi- 
tational interactions with a companion. There is strong 
observational support for the tidal circularization sce- 
nario. First, the current observations exhibit different 
orbital distributi ons for apparently single- and multiple- 
planet systems (jWright at al.ll2009f) . The so-called "3- 
day pileup" and a jump in planetary abundance beyond 
1 AU are not seen among multi-planet systems. It is dif- 
ficult to explain such a difference as a result of planet 
migration, since migration in the disk without strong 
planet-planet interactions is expected to lead to similar 
orbital distributions for both single and multiple planet 
systems. On the other hand, strong gravitational inter- 
actions between planets and/or stellar companions can 
disrupt the multiple-planet systems by ejecting planets 
or scattering them far away from one another. Thus, 
from such a scenario, it's expected that systems which 



have gone through violent dynamical interactions have 
only one planet close to the central star, while systems 
which did not experience such an event retain multiple 
planets close-by. Second, the observed distribution of 
stellar obliquities matches very well with the expecta- 
tion from Kozai migration (jTriaud et al.l[2010r i. Third, 
the observed inner edge of the orbital distribution is not 
at the Roche limit aR, which would be expected from 
planet migration, but rather at 2aR, which is naturally 
explained by circularization of a highly eccentric orbit 
while conserving total angular momentum (jFord fc Rasiol 
[2f)Ml . 

In this paper, we have further explored the possibility 
of forming close-in planets via tidal circularization of an 
eccentric orbit by studying the past and future evolution 
of eccentric transiting planets. We considered a broad 
range of tidal quality factors, ^ 10** — 10^, consistent with 
our limited theoretical understanding of tidal dissipation 
in both stars and planets. We have used a tidal model 
where Q' ex l/(crAi) as the simplest representation be- 
cause of the unknown character of tidal dissipation in 
stars and gaseous planets. Our choice is consistent with 
the constant time lag model ()Hutlll98l|) that is derived as 
a quadrupolar approximation of the tidal potential. Al- 
though different tidal models can change the time scales 
of the various processes, the qualitative properties of the 
evolution should be the same as other models. Another 
caveat here is that, for these close-in planets, the ratio 
R^,/a is relatively large and often ~ 0.1 (e.g., see Fig.IHl). 
Thus, the higher-order terms in the tidal potential can 
be important, as they are for M ars' satellites Phobos and 
Deimos (e.g. jBills etralll2005f) . 

In Section 2 we investigated the effects of (known or un- 
known) planetary/stellar companions on the evolution of 
observed close-in planets. Comparing secular timescales 
with GR precession timescales, we showed that the cur- 
rent and future evolution of most close-in planets is un- 
likely to be strongly affected by such companions. In 
Section 4.1 we re-examined the tidal stability of transit- 
ing systems, and confirmed that the majority of close- 
in planets are Darwin-unstable. The exceptions include 
CoRoT-3, CoRoT-6, HD 80606, and WASP-7 which are 
Darwin-stable within the current observational uncer- 
tainties. We also found that borderline cases like HAT- 
P-2 and WASP-10 are most likely Darwin-unstable, but 
could be Darwin-stable within the uncertainties. 

For clearly Darwin unstable systems (with Ltot ^ 
Lcrit), there are two possible evolutionary paths. When 
the tidal dissipation in the star dominates the evolution, 
all parameters but the planetary spin evolve on a similar 
timescale [t^ ^ Tq ^ r^^ ^ t^^ ) . On the other hand, 
when the tidal dissipation in the planet is non-negligible, 
the orbit tends to get circularized before the planet spi- 
rals all the way to the Roche limit (rg < Tq ^ r^, ^ Tui, ) . 
Although it is nontrivial to determine which evolutionary 
path each system will take without knowing the efficien- 
cies of tidal dissipation in the star and in the planet, there 
are some indications that dissipation in the planet dom- 
inates (rg < Ta) for most systems. First, there is clear 
evidence of eccentricity damping within ~ 0.1 AU. Most 
observed close-in planets are on nearly circular orbits, 
which are well-described by the traditional exponential 
eccentricity damping approximation. Second, our results 
in Section 5.2 suggest that we need to assume t^ < Ta 
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for most systems in order to explain the current obliq- 
uity distribution through Kozai migration. With stellar- 
dissipation dominated cases (re ~ Ta), we expect the 
distribution of current stellar obliquity would be signif- 
icantly different from the one expected from Kozai mi- 
gration. 

In Section 4.3, we showed that the lifetime of the plan- 
etary systems is largely determined by tidal dissipation 
in the star, while the circularization time is largely deter- 
mined by dissipati on in the pla,net. A lso, we confirmed 
the results by Barker fc Qgilvig l|20091 and showed that 
magnetic braking does not have a large effect for the fu- 
ture evolution except for systems with a rapidly rotating 
star. The minimum stellar tidal quality factor that al- 
lows a planet to survive for a certain age is similar for 
the cases with and without magnetic braking (see Fig. [7] 
and HI). 

In Section 5, we found that, generally speaking, the 
evolution history in the stellar-dissipation dominated 
tidal evolution is consistent with that expected from the 
planet migration origin of the close-in planets, with little 
change in semi-major axis and eccentricity over the stel- 
lar age. On the other hand, the evolution history in the 
planetary-dissipation dominated evolution is consistent 
with that expected from scattering/Kozai-cycle origin, 
with initial orbits being wide and eccentric. The latter 
case agrees with the results in I Jackson et al.l |2008). We 
also showed that, when the effects of tidal dissipation 
in the star is comparable to, or larger than that in the 



planet, the stellar obliquities could be damped from high 
> 40° to low ^ 2° values within the currently observed 
stellar ages (see Fig. [151 and [TO)) . 

Overall, our results for the tidal evolution of eccen- 
tric transiting planets is consistent with the formation 
path that involves circularization of an initially eccen- 
tric orbit. The distribution of system parameters seems 
to imply that this mechanism dominates for close, single 
planet systems, but multiplanet systems are more consis- 
tent with the migration to account for the close members. 



This work was supported by an Astronomy Center for 
Theory and Computation Prize Fellowship (to S. M.) at 
the University of Maryland, and by NSF Grant AST- 
0507727 (to F. A. R.) at Northwestern University. Par- 
tial support was also provided by the Center for Inter- 
disciplinary Exploration and Research in Astrophysics 
(CIERA) at Northwestern University. S. J. Peale is sup- 
ported in part by NASA under grant NNG06GF42G 
in the Origins of Solar Systems Program. S. M. also 
thanks the support from the Kavli Institute for Theoret- 
ical Physics at UC Santa Barbara while part of this work 
was done. We thank Benjamin Levrard, Giles Chabrier, 
Brian Jackson, Amaury Triaud, Dong Lai, Yoram Lith- 
wick, Smadar Naoz for useful discussions, and an anony- 
mous referee for valuable comments on the manuscript. 
S. M. thanks Genya Takeda for making the time com- 
parison code used in Fig. [1] and [2] available. 



REFERENCES 



Adams, F. C. & Laughlin, G. 2003, Icarus, 163, 290 

— . 2006, ApJ, 649, 992 

Alexander, M. E. 1973, Ap&SS, 23, 459 

Alonso, R., Auvergne, M., Baglin, A., OUivier, M., Moutou, C, 
Rouan, D., Deeg, H. J., Aigrain, S., Almenara, J. M., Barbieri, 
M., Barge, P., Benz, W., Borde, P., Bouchy, F., de La Reza, 
R., Deleuil, M., Dvorak, R., Erikson, A., Fridlund, M., Gillon, 
M., Gondoin, P., Guillot, T., Hatzes, A., Hebrard, G., Kabath, 
P., Jorda, L., Lammer, H., Leger, A., Llebaria, A., Loeillet, B., 
Magain, P., Mayor, M., Mazeh, T., Piitzold, M., Pepe, F., Pont, 
F., Queloz, D., Rauer, H., Shporer, A., Schneider, J., Stecklum, 
B., Udry, S., & Wuchterl, G. 2008, A&A, 482, L21 

Barge, P., Baglin, A., Auvergne, M., Rauer, H., Leger, A., 
Schneider, J., Pont, F., Aigrain, S., Almenara, J., Alonso, R., 
Barbieri, M., Borde, P., Bouchy, F., Deeg, H. J., La Reza, D., 
Deleuil, M., Dvorak, R., Erikson, A., Fridlund, M., Gillon, M., 
Gondoin, P., Guillot, T., Hatzes, A., Hebrard, G., Jorda, L., 
Kabath, P., Lammer, H., Llebaria, A., Loeillet, B., Magain, P., 
Mazeh, T., Moutou, C., OUivier, M., Patzold, M., Queloz, D., 
Rouan, D., Shporer, A., & Wuchterl, G. 2008, A&A, 482, L17 

Barker, A. J. & Ogilvie, G. L 2009, MNRAS, 395, 2268 

Bills, B. G., Neumann, G. A., Smith, D. E., & Zuber, M. T. 2005, 
Journal of Geophysical Research (Planets), 110, 7004 

Bodenheimer, P., Lin, D. N. C., & MardUng, R. A. 2001, ApJ, 548, 
466 

Chandrasekhar, S. 1987, Ellipsoidal figures of equilibrium (New 
York : Dover, 1987.) 

Ghatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, 
ApJ, 686, 580 

Darwin, G. H. 1879, proc. Roy. Soc. 

Dobbs-Dixon, L, Lin, D. N. G., & MardUng, R. A. 2004, ApJ, 610, 
464 

Ecuvillon, A., Israelian, G., Santos, N. G., Mayor, M., & Gilli, G. 
2006, A&A, 449, 809 

Eggenberger, A., Udry, S., & Mayor, M. 2004, A&A, 417, 353 

Fabrycky, D. & Tremaine, S. 2007, ApJ, 669, 1298 

Fabrycky, D. C. & Winn, J. N. 2009, ApJ, 696, 1230 

Ferraz-Mello, S., Rodriguez, A., & Hussmann, H. 2008, Gelestial 
Mechanics and Dynamical Astronomy, 101, 171 

Ford, E. B., Kozinsky, B., & Rasio, F. A. 2000, ApJ, 535, 385 



Ford, E. B. & Rasio, F. A. 2006, ApJL, 638, L45 

Gaudi, B. S. & Winn, J. N. 2007, ApJ, 655, 550 

Gillon, M., Smalley, B., Hebb, L., Anderson, D. R., Triaud, 

A. H. M. J., Hellier, C., Maxted, P. F. L., Queloz, D., & Wilson, 

D. M. 2009, A&A, 496, 259 
Goldreich, P. 1963, MNRAS, 126, 257 
Goldreich, P. & Tremaine, S. 1980, ApJ, 241, 425 
Gonzalez, G. 1997, MNRAS, 285, 403 

Gu, P., Lin, D. N. G., & Bodenheimer, P. H. 2003, ApJ, 588, 509 
Hale, A. 1994, AJ, 107, 306 
Hirano, T., Narita, N., Shporer, A., Sato, B., Aoki, W., & Tamura, 

M. 2010, ArXiv e-prints 
Holman, M., Touma, J., & Tremaine, S. 1997, Nature, 386, 254 
Hut, P. 1980, A&A, 92, 167 
— . 1981, A&A, 99, 126 

Ibgui, L. & Burrows, A. 2009, ApJ, 700, 1921 
Jackson, B., Barnes, R., & Greenberg, R. 2009, ArXiv e-prints 
Jackson, B., Greenberg, R., & Barnes, R. 2008, ApJ, 678, 1396 
Johns-KruU, G. M., McCuUough, P. R., Burke, C. J., Valenti, J. A., 

Janes, K. A., Heasley, J. N., Prato, L., Bissinger, R., Fleenor, 

M., Foote, C. N., Garcia-Melendo, E., Gary, B. L., HoweU, P. J., 

MaUia, F., Masi, G., & Vanmunster, T. 2008, ApJ, 677, 657 
Johnson, J. A., Winn, J. N., Albrecht, S., Howard, A. W., Marcy, 

G. W., & Gazak, J. Z. 2009, PASP, 121, 1104 
Kiseleva, L. G., Eggleton, P. P., & Mikkola, S. 1998, MNRAS, 300, 

292 
Kozai, Y. 1962, AJ, 67, 591 

Lai, D., Rasio, F. A., & Shapiro, S. L. 1994, ApJ, 423, 344 
Lainey, V., Arlot, J., Karatekin, O., & van Hoolst, T. 2009, Nature, 

459, 957 
Lammer, H., Odert, P., Leitzinger, M., Khodachenko, M. L., 

Panchenko, M., Kulikov, Y. N., Zhang, T. L., Lichtenegger, 

H. I. M., Erkaev, N. V., Wuchterl, G., Micela, G., Penz, T., 

Biernat, H. K., WeingriU, J., Steller, M., Ottacher, H., Hasiba, 

J., & Hanslmeier, A. 2009, A&A, 506, 399 
Lammer, H., Selsis, F., Ribas, I., Guinan, E. F., Bauer, S. J., & 

Weiss, W. W. 2003, ApJL, 598, L121 
LecaveUer Des Etangs, A. 2007, A&A, 461, 1185 
Leconte, J., Chabrier, G., BarafFe, L, & Levrard, B. 2010, A&A, 

516, A64+ 



19 



Levrard, B., Correia, A. C. M., Chabrier, G., BaraflFe, I., Selsis, F., 

& Laskar, J. 2007, A&A, 462, L5 
Levrard, B., Winisdoerffer, C, & Chabrier, G. 2009, ApJL, 692, 

L9 
Lidov, M. L. 1962, P&SS, 9, 719 
Love, A. E. H. 1944, A Treatise on the Mathematical Theory of 

Elasticity (Mineola, NY: Dover) 
Mardling, R. A. & Lin, D. N. C. 2004, ApJ, 614, 955 
Matsumura, S., Takeda, G., & Rasio, F. A. 2008, ApjL, 686, L29 
Matsumura, S., Thommes, E. W., Chatterjee, S., & Rasio, F. A. 

2010, ApJ, 714, 194 
McLaughlin, D. B. 1924, ApJ, 60, 22 

Miller, N., Fortney, J. J., & Jackson, B. 2009, ApJ, 702, 1413 
Murray, N., Hansen, B., Holman, M., & Treinaine, S. 1998, Science, 

279, 69 
Nagasawa, M., Ida, S., & Bessho, T. 2008, ApJ, 678, 498 
Ogilvie, G. L & Lin, D. N. C. 2004, ApJ, 610, 477 
— . 2007, ApJ, 661, 1180 

Ohta, Y., Taruya, A., & Suto, Y. 2005, ApJ, 622, 1118 
Paczynski, B. 1971, ARA&A, 9, 183 
Pal, A., Bakos, G. A., Torres, G., Noyes, R. W., Fischer, D. A., 

Johnson, J. A., Henry, G. W., Butler, R. P., Marcy, G. W., 

Howard, A. W., Sipocz, B., Latham, D. W., & Esquerdo, G. A. 

2010, MNRAS, 401, 2665 
Papaloizou, J. C. B. & Larwood, J. D. 2000, MNRAS, 315, 823 
Peale, S. J. 1986, Orbital resonances, unusual configurations and 

exotic rotation states among planetary satellites (Burns, J. A. & 

Matthews, M. S., U. Ariz. Press, Tucson,m AZ), 159-223 
Pont, F., Endl, M., Cochran, W. D., Barnes, S. L, Sneden, C, 

MacQueen, P. J., Moutou, C, Aigrain, S., Alonso, R., Baglin, 

A., Bouchy, F., Deleuil, M., Fridlund, M., Hebrard, G., Hatzes, 

A., Mazeh, T., & Shporer, A. 2010, MNRAS, 402, LI 
Rasio, F. A., Tout, C. A., Lubow, S. H., & Livio, M. 1996, ApJ, 

470, 1187 
Rossiter, R. A. 1924, ApJ, 60, 15 
Skumanich, A. 1972, ApJ, 171, 565 
Sterne, T. E. 1939, MNRAS, 99, 451 



Takeda, G., Kita, R., & Rasio, F. A. 2008, ApJ, 683, 1063 
Torres, G., Winn, J. N., & Holman, M. J. 2008, ApJ, 677, 1324 
Triaud, A. H. M. J., Collier Cameron, A., Queloz, D., Anderson, 

D. R., Gillon, M., Hebb, L., Hellier, C, Loeillet, B., Maxted, 
P. F., Mayor, M., Pepe, F., PoUacco, D., Segransan, D., Smalley, 
B., Udry, S., West, R. G., & Wheatley, P. J. 2010, ArXiv e-prints 

Trilling, D. E., Benz, W., Guillot, T., Lunine, J. L, Hubbard, W. B., 

& Burrows, A. 1998, ApJ, 500, 428 
Ward, W. R. 1997, Icarus, 126, 261 
Winn, J. N., Fabrycky, D., Albrecht, S., & Johnson, J. A. 2010a, 

ApJL, 718, L145 
Winn, J. N., Howard, A. W., Johnson, J. A., Marcy, G. W., 

Gazak, J. Z., Starkey, D., Ford, E. B., Colon, K. D., Reyes, F., 

Nortmann, L., Dreizler, S., Odewahn, S., Welsh, W. F., Kadakia, 

S., Vanderbei, R. J., Adams, E. R., Lockhart, M., Crossfield, I. J., 

Valenti, J. A., Dantowitz, R., & Carter, J. A. 2009a, ApJ, 703, 

2091 
Winn, J. N., Johnson, J. A., Albrecht, S., Howard, A. W., Marcy, 

G. W., Crossfield, I. J., & Holman, M. J. 2009b, ApJL, 703, L99 
Winn, J. N., Johnson, J. A., Fabrycky, D., Howard, A. W., Marcy, 

G. W., Narita, N., Crossfield, I. J., Suto, Y., Turner, E. L., 

Esquerdo, G., & Holman, M. J. 2009c, ApJ, 700, 302 
Winn, J. N., Johnson, J. A., Howard, A. W., Marcy, G. W., Bakos, 

G. A., Hartman, J., Torres, G., Albrecht, S., & Narita, N. 2010b, 

ApJ, 718, 575 
Winn, J. N., Johnson, J. A., Howard, A. W., Marcy, G. W., 

Isaacson, H., Shporer, A., Bakos, G. A., Hartman, J. D., & 

Albrecht, S. 2010c, ArXiv e-prints 
Wright, J. T., Upadhyay, S., Marcy, G. W., Fischer, D. A., Ford, 

E. B., & Johnson, J. A. 2009, ApJ, 693, 1084 
Wu, Y. 2005, ApJ, 635, 688 

Wu, Y. & Murray, N. 2003, ApJ, 589, 605 

Wu, Y., Murray, N. W., & Ramsahai, J. M. 2007, ApJ, 670, 820 

Zhang, K. & Hamilton, D. P. 2008, Icarus, 193, 267 

Zhou, J.-L. & Sun, Y.-S. 2003, ApJ, 598, 1290 



Table 1. Data are taken from 'http : //exoplanet . eu/ ', 



to 

o 



Planet Name 



Mp 
[Mj] 



[Rj] 



[AU] 



[Mq] 



[Re] 



II sin « 
[km/s] 



A 
[degrees] 



Age 
[Gyr] 



References 



CoRoT-1 b (GOV) 
CoRoT-2 b (G7V) 
CoRoT-3 b (F3V) 
CoRoT-4 b (FOV) 
CoRoT-5 b (F9V) 
CoRoT-6 b (F9V) 
CoRoT-7 b (G9V) 
CoRoT-7 c (G9V) 
GJ 1214 b (M) 
GJ 436 b (M2.5) 
HAT-P-1 b (GOV) 
HAT-P-11 b (K4) 
HAT- P- 12 b (K4) 
HAT- P- 13 b (G4) 
HAT-P-13 c (G4) 
HAT-P-2 b (F8) 
HAT-P-3 b (K) 
HAT-P-4 b (F) 
HAT-P-5 b (G) 
HAT-P-6 b (F8) 
HAT-P-7 b (F6V) 
HAT-P-8 b (F) 
HAT-P-9 b (F) 
HD 149026 b (GOIV) 
HD 17156 b (GO) 
HD 189733 b (Kl-2) 
HD 209458 b (GOV) 
HD 80606 b (G5) 
Kepler-4 b (GO) 
Kepler-5 b (?) 
Kepler-6 b (F) 
Kepler-7 b (F-G) 
Kepler-8 b (F8IV) 
Lupus-TR-3 (KIV) 
OGLE-TR-10 (G/K) 
OGLE-TR-111 (G/K) 
OGLE-TR-113 (K) 
OGLE-TR-132 (F) 
OGLE-TR-182 (G) 
OGLE-TR-211 (F7-8) 
OGLE-TR-56 (G) 
TrES-1 (KOV) 
TrES-2 (GOV) 
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Oi oq-t-0.S2 

Zl.ZO_Q gg 
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9 Qfj+0.34 

n 0151+00025 

U.Uiai_o 0025 

0264+00028 

U.U^04_o 0028 
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0720+00025 

U.U/za_o 0025 

n rrqn+O.OSO 

nnsi+0009 

U.0»l_o.oo9 

211+0012 
U.^il_0.oi2 
n = 0+0.029 
U.(500_o 046 

15.21;:° 

9-091^:^4 

596+0024 
U.OyO_o 026 

n fis+0-04 
U.b8_o.o4 

1 06+0-" 
l-'-'O-O.ll 

1 057+01" 
1 776+0077 
^•' '"-0.049 
1 C9+0.18 

j+0.09 

'-0.09 

n Q^n+0.022 

u.ooy_o 021 

q 99 + 0.08 

■^•^^-0.08 
1 iqs+0-022 
1.100_o 025 

0.685«:°;i 

077+0.012 
'^•'-" '-0.012 
9 -, ,.+0.056 
^•^^^-0.059 
669+0025 

U.DD9_o 030 
433+0.040 
u.y:00_o 041 

o.603t;;;^ 

n 81+0.18 

62+0-14 
'^•°^-0.14 

55+0-10 
'J-^^-o.io 

1 26+016 
l-^"-0.16 

O+0.14 

'-0.13 

101+0-15 
1-'J1-0.15 

q+0.20 

^-0.20 
1 30+0.18 

752+0-047 
U.(OZ_o 046 

1.200+°-°5l 



0.78^ 



4.20^ 



1.18'* 



1.03^ 



1 40+0.08 
l-*y-0.08 
1 465+0029 
i.4DO_o 029 

n Q034+0-05* 
U.MM01_o 058 
1 lon+ooe 
i.iyu_o 05 

1 qos+0.046 
l.ooo_o,o47 



1.166] 



n+0.035 
-0.035 
n. 1 c:+0.008 
'J-^3-0.008 



0.239 

n q7f;7+0.0082 
U.O(D(_o,oo92 



1.242"' 



1+0.053 

^=-0.053 

n 422+0-014 
U.4ZZ_o 014 

n 050+0.029 
u.Moy_o 021 

1.281 



+0.079 



-0.053 



1 1=7+0.073 
-'-•-'-'J' -0.062 

0.899t°-.°« 

1 97+0.05 
l-^'-0.05 

1 26+0-05 
1-^0-0.05 
1 qqn+0.061 
l-oou_o 061 
1 363+0-195 
i.ODO_o 087 

1 50+0-08 
l-OU_o 06 

1 40+0-06 

i-^'J-o.oe 

654+0-060 
U.004_o 045 

1 02+0-08 
^^•'-'^-0.08 

1 170+0. 016 
l-l'O-0.023 
1 350+0016 
l.OO3_0 019 

074+0.030 
u.»('±_o 030 
p, q^„+0.019 
U.OO/_o 019 

1 431+0041 
l.'iOl_o 052 
1 q9q+0.026 
i.OZO_o 029 
1 470+0.050 
l-478_oo5i 

i-4i9«:°f8 

0.89tH? 

1 9 =+0.14 
1-^0-0.12 
1 051+0-057 
i.U01_o 052 

1 093+0-028 

l.U30_o 019 

1 20+0-15 

1 ,q+0.24 
^•4'J-0.08 

1 36+0-18 
l-JO-0.09 
1 qcq+0.092 
i.ouo_o 090 
1 067+0-022 

i.UD/_o 021 

1 224+0-041 
l.zzq:_o 041 



0.0254"* 



0.09: 
0.04947"^ 



.,+0.0004 
-0.0004 

n nosi +0-0009 

O.U2»l_ooo09 
05694+0-00096 
U.UaDM4_o 00079 
3+0.001 
-0.001 

f-0. 00026 
-0.00029 

0.0855«:°;^« 

n ni79+000029 
U.Ul(Z_o 00029 

0.046 
0.0143 

n n9S79+000029 
U.UZO/z_o 00026 
0553+0-0012 
U.UOOO_o 0013 
0530+0-0002 

U.UOOU_o 0008 
O O'^«/l+00003 
O.U384_ooo03 

0427+0-0006 
U.U4Z/_0 0012 

1 1OS+0-018 
l.i»»_0.033 

n Ofi«7S+0-00068 
0-UD»/»_o. 00068 
n nq«ao+0. 00060 
O.U38»2_o 00077 

0.0446t°:°?^^ 

-+0.00076 
^-0.00076 
05235+0-00087 
U.UOZOO_o 00087 
n 0^77+00005 
U.UJ^ /_0 0005 
O O/1S7+0-0026 
O.U48/'_o.oo26 
n 0=3+0.002 
u.uoo_o 002 
04313+0-00065 
U.U4cilci_o 00056 
1614+0-0022 
U.10i4_o 0022 
03120+0-00027 
U.UOiZU_o 00037 
04707+0-00046 
U.Ut(U(_o 00047 
4614+0-0047 
U.1Diy:_o 0047 
0456+0-0009 
U.uy:OD_o 0009 
05064+0-00070 
U.UOUO'l_o 00070 
04567+0-00055 



0.04075^ 



-0.00046 
0.06224lHgr4 
0.0483t°:^^ 
0464+00007 
U.U^Dy:_o 0007 

0.0434t°:°«;^ 
0.04689j:°:~j^7 

O O99SO+0-00016 
0.02289_o 00015 
03035+0-00057 
U.UOUOO_o 00053 

051+0-001 
u.uai_o 001 

051+0-001 

U.U01_o 001 
n n9ctOQ+0. 00046 
O.U2ci»d_o 00051 
03925+0-00056 
u.uo»Zd_o 00060 
O 0^558+0-00070 
U.U3ODO_0 00077 



+0.1 



0.015 
0.005 



0.1 
+0.09 
0.04 



(fixed) 

(fixed) 

0.008 



0.09 

<0.1 





< 0.27 

14+0-01 
U.14_0.oi 

< 0.067 

O 1OS+0-046 
0.198_o.o46 







,+0.009 



0.021_o 009 
691+0-018 

u.oyi_o 018 

n =17-1+0.0033 
U.01(l_o 0033 







(fixed) 

(fixed) 

(fixed) 

(fixed) 



O 6SOI +0-0019 

0.bSUi_o.ooi9 

0041+0-0025 

U.UU41_0 0020 



0.932861 

(fixed) 

< 0.024 

(fixed) 

(fixed) 

(fixed) 

(fixed) 





















0.00055 
0.00055 



0.95 

0.97 

1.359 

1.16 

1.00 

1.05 

0.93 



+0.15 

0.15 
+0.06 

0.06 

+0.059 
0.043 
+0.03 

0.02 
+0.02 

0.02 
+0.05 

0.05 
+0.03 

0.03 



0.157 

0.452 

1.133 

0.81 

0.733 

1.22 



+0.019 
0.019 

+0.014 
0.012 

+0.075 
0.079 
+0.02 

0.03 

+0.018 

-0.018 
+0.05 

0.10 



1 36+0-04 
l-JD-0.04 

O9S+0-044 
0.928_o.o54 

1 26+0-06 
l-^"-0.14 

1 160+0-062 
1.10U_o 062 
1 9Q+0.06 
1-^^-0.06 
-, 47+0.08 
l-^'-0.05 
-I 90+0.04 
l-2»_o.o4 
1 90+0.13 
l-^*'-0.13 
1 9Q-+0.060 
l-^y^-0.050 
1 94+0-03 
1-^^-0.03 
n O9q+0.022 
U.OZO_o 029 
1 110+0.033 
l.llj_0 033 
1 05+0.032 
1-'J3-0.032 
1 99q+0.053 
l.zzo_o 091 
1 q7.+0.040 
^-•^'^-0.059 

1 209+0-044 
l.zu»_o 038 
1 347+0.072 
l.oy:(_o 054 
1 9iq+0.067 
l.zio_o 063 

n 07+0.04 

0-»'_o.04 
1 14+0.10 
l--^*-0.12 

0=9+0.058 
U.OOZ_o 052 
Q 770+0.017 
"■' '^-0.015 

1 305+0-075 
1.0UO_o.067 
1 14+0.05 
1-1^-0.05 

-, qq+0.05 

l-'->'J-0.05 
-, 990+0.072 
l.zzo_o 078 
n S7S+0-038 
0-»'»_0.040 
O oS'^+0-059 
u.yoo_0 og3 



1 -,1+0.05 
1-11-0.05 

009+0-018 

u.yuz_o 018 

1 540+0-083 
1.04U_o 078 
-, -17+0.01 
l-l'-0.03 

1 1S6+0-04 
l.i8b_o.o4 
1 095+0.026 
l.UZO_o 026 

n 07+0.04 

0.»/_o.04 



0.211 

0.464 

1.135 

0.75 

0.701 

1.56 



+0.0097 
0.0097 

+0.009 
0.011 

+0.048 
0.048 
+0.02 

0.02 

+0.02 

-0.01 
+0.08 

0.08 



1 f;4+0.09 
l-°*-0.08 
n oqq+0.034 
U.OOO_o 044 
1 =0+0.07 

i-oy-0.07 

1 167+0-049 
l.iO(_0 049 

1.8410-23 

1 =0+0.08 
l-0o_o 06 
-, q9+0.07 
-'-■'J^-0.07 

1.368_o 083 
-1 44+0-08 
l-^*-0.08 
766+0-007 

U.(OD_o 013 

-I 1 ==+0.014 
1.100_o 016 

O6S+0-028 
0-9b8_oo28 

1 .07+0.071 
l.'iO(_o 084 
1 7Qq+0.043 
l.(yo_o 062 
1 301+0.017 
l.oyi_o 034 
1 o/iq+0.048 
l-o'±o_o 066 
1 4Sfi+0-053 
1.48b_oo62 

n 09+0.05 

^■*'^-0.05 

-, 17+0.13 

l-l'-O.ll 
n sqi +0-045 
O.»di_oo40 

771+0.020 

"■' '^-0.011 

1 q9+0.17 
-^-'-'^-0.12 
1 14+0.23 
1-1^-0.06 

1 64+0-21 
l-°*-0.07 

1 qf,q+0.089 

l.ooo_o 086 

0.8071H;^ 



5.2 
11.85 
17.0 
6.4 



+1.0 
1.0 

+0.50 
0.50 
+ 1.0 
1.0 
+ 1.0 
1.0 



7.6 



+1.0 
1.0 



< 2.0 

0.52 

3.75 

1.5 

0.5 

2.9 



+0.05 
0.05 

+0.58 
0.58 
+ 1.5 

1.5 
+0.4 

0.4 
+ 1.0 

1.0 



20.8 



+0.3 



-0.3 

p;+0.5 
-0.5 
c;+0.5 
-0.5 
K+1-5 
-1.5 
7+1.0 
-1.0 
i+0.5 
-0.5 
-+0.5 
-0.5 
3+1-0 
-1.0 
1+2.1 
-0.6 
a+0.31 
-0.31 
.,+0.017 
-0.067 
1+0.16 
-0.16 
1+0.44 
-0.22 
1+1.0 
-1.0 
a+1-0 
-1.0 
-1.0 
-1.0 
,+0.5 
-0.5 

10.5+0-^ 



0.5: 
5.5J 

2.6;; 

8.7I 

3.8^ 

11.5; 

11. 9J 

6.2^ 

4.18! 

3.316! 

4.70J 

1.12J 

2.2"* 

4.) 

3.0! 

4.2"' 



-0.7 



,+2 
^-2 



3 

1.3 

1.0 



+0.3 
0.3 

+0.6 
0.6 



77+11 
"-11 




BargeOS, PontlO 


7 2+4-5 
'■^-4.5 


-0.2-4 


AlonsoOS, Bouchy08 


37 R+22.3 
— d'-D-lO.O 


1.6-2.8 


Triaud09, Deleuil08 




1+1.0 
1-0.3 


MoutouOS 




5.5-8.3 


RauerOg 




2.5-3.3 


FridlundlO 




1.2-2.3 


QuelozOO 
Queloz09 




3-10 


Charbonneau09 




6.01^:1; 


TWC08, LWC09 


3 7+2.1 
■^-'-2.1 


9 7+2.5 
^■'-2.0 


TWC08, Johnson08 




6.5lt:? 


BackoslO 




2.5l|g 


Hartman09 




5.01^:5 


Bakos09a 
Bakos09a 


1 9+13.4 
1-^-13.4 


2.61SJ 


PallO, Winn07 




l-5l?1 


TWC08, Torres07 


182.51^:1 


4.2l§:^ 

2.6i;:^ 

9 q+0.5 

9 9 + 1.0 


Kovacs07 
Bakos07 
Noyes08 
Pal08, Winn09c 




3.4i;:0 


Latham09 


-i2i;5 
10.015-; 

o =+0.32 
-U.»tl_o28 

-4.41;-.; 

531^2^ 


i.6i;i 

i-9lS:^ 

3.061°;?* 

q 1+0.8 
■^■1-0.7 

i.6i;:? 

4.5i;-5 


Shporer09 

TWC08, Wolf07 

Barbieri09, Narita09, Winn09a 

Triaud09, TWC08 

TWC08, Winn05 

Winn09b 

BoruckilO 


-26.9l*:^ 


30+0-6 

q O+l'O 

q = + 1.0 

3.84i;-5 


KochlO 

DunhamlO 

LathamlO 

JenkinslO 

WeldrakeOS 




q 9 + 4.0 

0+5.2 

»-»-6.6 

13.2+0-8 

1 2+1-5 
1-^-1.1 


TWC08, KonackiOS 

TWC08 

TWC08 

TWC08 

Pont08 

UdalskiOS 


301-1 

-9.oi;-° 


q 9 + 1.0 
q 7+3.4 

5-ol^:I 


TWC08, Konacki03 
TWC08, Narita07 
TWC08, WinnOS 
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Mp 


Rp 


a 


e 


M* 


R* 


vsini 


A 


Age 


References 




[Mj] 


\Rj] 


[AU] 




[Mq] 


[Re] 


[km/s] 


[degrees] 


[Gyr] 




TrES-3 (G) 


1 QQo+0.062 

l.aOO_g gg3 


, qi 9+0.033 


0.02272J:g;~S^ 





OQI ^+0-021 
U.3J-0_o 031 


019+0.014 

U-»i^_0.025 


1 t;+i-0 




0.61^4 


TWC08 


TrES-4 (F) 


0.920+0-°^3 


1 yc;-, +0.064 
^^•'0^-0.062 


0.05092+_'oZll 





1 OQ4+0.060 
4-oy4_o 056 


1 C.1C+0.065 
1.8ib_o.o62 


c: + 1.2 
».ti_l 2 


6.3lt^7 


2.9ir4 


TWC08, NaritalO 


WASP-1 b {F7V) 


U.aiO_Q ggg 


1 CI 4+0.052 


0.03957lg:S 





1 ■HOI +0-049 
l.dUi_o.o47 


1 r,i 7+0.052 
l--Jl-'_0.045 


c 7q+0.35 
0-'y_0.35 




0+0.6 


TWC08, Stempels07 


WASP-10 b (K5) 


2.96l«:f, 


1 90+0. 077 


0.0369t°:g°l^ 


0.0591HJ1 


f. 7f,q+0.068 
U.CUcl_o.080 


77^+0.043 
U.MO_o 040 


<6 




n s+O-2 

0-8-0.2 


Christian09 


WASP- 11 b (K3V) 


oAsrtVoll 


1 nn^+o''32 


nA-i^+ooooe 

U.U4cit)_o 0006 


(fixed) 


n OO+0.03 
^■°'^-0.03 


7q+0.02 
^■'^-0.02 


K+0.2 
U.&_o.2 




7-9tli 


Bakos09b 


WASP-12 b (GO) 


i-4ilS:iS 


, 7Q+0.09 
1- '9-0.09 


099Q+00008 
U.U^^9_ooo08 


0.049lHi^ 


1 qc:+0-14 
-■-■■^3-0.14 


-1 c:7+0.07 


< 2.2II j 




2tl 


Hebb09 


WASP-13 b (GIV) 


0.461H! 


1 91+0.14 


0.0527^0;°°}^ 


(fixed) 


1 03+0-" 
J^-UJ_o.o9 


-I 04+0.13 


<4.9 




8.5tll 


Skillen09 


WASP-14 b (F5V) 


7 q --,+0.508 
' -^^^-0.496 


1 901+0.075 
i.ZSi_o 082 


0.036tHgl 


oqi+0003 
U.U9i_o.o03 


-, 911+0.127 
^■^^J^-0.122 


1 QOfi+0-066 
i.30b_o.o73 


4Q+1-0 


-33.lt?i 


-0.5-1.0 


Joshi09, Johnson09 


WASP-15 b (F5) 


o.542«;°i° 


1 /19S+0-077 
l-4iS_o.o77 


0.0499t0:«gl| 


(fixed) 


1 1S+0-12 
^■^»-0.12 


1.477tH?^ 


^■27tlil 


-139.6t^-J 


3-9l?:3 


West09, TriaudlO 


WASP-16 b (G3V) 


0.855t°;°« 


1 nns+0-083 

l.UU8_oo60 


0.0421«:rig 


(fixed) 


-, ^199+0.074 
l.U^2_o 129 


0.9461HS 


3.0+Jd 




2-3+2'2 


Lister09 


WASP- 17 b (F6) 


OAmtVoll 


1 74+0.26 
^•'4-0.23 


0.0501«:ri^ 


0.129l«:S°| 


1 90+0-12 
l-^'J-0.12 


1 00+0.20 

l--'3»_o.l8 


W.14t°il 


-147.3t^-J 


3.oto;^ 


AndersonlO, TriaudlO 


WASP-18 b (F9) 


^0.30t°oil 


1 ip,f,+0.072 


0.02026lg:rS 


0.0092^0.^21 


-, 9C+O.13 
l-^ti-0.13 


1 9lfi+0.067 
1--^1-O_0.054 


14.67+0-«i 


5.ot^3:? 


0.5-1.5 


HeUier09a, TriaudlO 


WASP-19 b (G8V) 


l-14tH? 


1 90+0.07 
i-^»_0.07 


01fi4+0"005 
U.UiD4_o 0006 


o.02«;«? 


qc+O-lO 
^■^3-0.10 


qq+0.05 
U.30_o 04 


4±^ 




^1 


HebblO 


WASP-2 b (KIV) 


0.915^2:^3 


1 f.7, +0.080 
4-U'i_o.o83 


n n-:!1'?S+000130 
U.U3138_o.ooi54 





0.89tri2 


n szin+0-062 

0.840_o.o65 


0.99+_lil 


-153+15 


5.6l^:^ 


TWC08, TriaudlO 


WASP-3 b (F7V) 


1 7f;+0.08 


1 01 +0.07 
i-Ji_0.14 


f.p,q, 7+0.0005 
U.Udi/^_0 0010 





i.24«-:;? 


1 q-i+0.06 
-'--'-'-'--0.12 


ISAtli 


15+10 


0.7-3.5 


Pollacco08, SimpsonOg 


WASP-4 b (G7V) 


1 r,i+0.13 

^•^^-0.08 


1 qo4+0-''^* 

i.JU4_oo42 


099c;c;+0. 00095 


(fixed) 


0-S0_o.o7 


n 070+0.036 

U.S/3_o.027 


9 -,4+0.38 


4+34 
4-43 


c 9+3.8 


Gillon09b, TriaudlO 


WASP-5 b (G4V) 


1 rro + 0.13 
i-t>»_0.10 


1.087tH?? 


n 09fi7+0"012 
U.UZD/^_0 0008 


0.038tH?^ 


qfi+O-13 

u.yo_o.o9 


1 099+0-056 
l.UZy_0.069 


91+0.34 
■^•^'^-0.35 


12.4t?-o 


5.4lt^ 


Gillon09b, TriaudlO 


WASP-6 b (G8V) 


t^oq+oois 


1 994+0.051 


0491+0-0008 
U.U42i_oooi3 


0.054lHi^ 


n ssn+O-050 

O.S»U_o.o80 


n 070+0.025 
0.870_o.o36 


i-4tl;g 


0.2010.^2 


ntl 


Gillon09a 


WASP-7 b (F5V) 


0.96lg:ii 


Oqic+0.046 

u.yit)_o 040 


00618-0. 0033 


(fixed) 


-I 90+0.09 


1 900+0.059 
1.2db_o.o46 


irtl 






HeUier09b 


XO-1 b (GIV) 


0.918«;g?^ 


1 9nc+0.047 
i.ZUO_o 042 


0.04928l°:»g99 





-, 097+0.057 

1-U^I'_0.061 


0.9341HS 


1 -,1+0.67 
1-11-0.67 




1 0+3-1 
1-0-0.9 


TWC08, McCuUoughOe 


XO-2 b (KOV) 


0.566lO;°f5 


0.983_oo28 


0.036840^00040^3 





074+0.032 

U.y/4_oo34 


0.97ltH^^ 


1 4+0.3 
1-4-0.3 




tr + 2.8 
t).»_2 3 


TWC08, Burke07 


XO-3 b (F5V) 


13.25«;^* 


1 Qc:+0.16 


n n/i'7C+0.0005 
0.0476_ooo05 


o.26otnr7 


1 4-1+0.03 
-"-■4-^-0.05 


9 10+0. 04 
^■l'^-0.05 


18.54«-1? 


37.3ti^-^ 


2.69«:1* 


Johns-KruUOS, Winn09d 


XO-4 b (F5V) 


-I 7r| + 0.20 

^^•'^-0.20 


1 qi+0.048 

J^-^^-0.048 


oc;c;c:+0.0011 
u.uooo_o 0011 


(fixed) 


1 09+0.02 
^■■'^-0.02 


1 ^fi+0-05 


0+0.5 

°-»-0.5 




9 -1+0.6 
^-1-0.6 


McCuUoughOS 


XO-5 b (G8V) 


1 O^Q+0''28 
i.UOy_Q 028 


1 inq+0.050 
i.iU9_oo50 


n nziss+0-0006 
0.0488_o.ooo6 





n SH+003 
0-8O-0.03 


1 nH+004 
l-U8_o.o4 


7+0.5 
0-'-0.5 




-,4 0+2.0 

14.8_2.o 


Pal09 



Note. — Column 1 — planet's name and the stellar spectral type inside the bracket, Column 2 — planetary mass, Column 3 - planetary radius, Column 4 - semi-major axis, Column 5 — eccentricity. 
Column 6 — stellar mass, Column 7 - stellar radius, Column 8 — projected stellar rotational velocity. C olumn 9 - projecte d stellar obli quity, Column 10 — stellar age. Column 11 — references, where 
we indicate the references by the first author name and the published year as, for example, BargeOS for lBarge et aU JSOOSl) . TWH08 is lTorres et al.| l |2008D 
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Table 2. Critical conditions for tidal instability, as well as the 
Roche limit for each system in Table [T] 



Planet Name 


n/uif 


n/(aj, cose) 


Prot,* [days] 


Ltot/ Lc 


-'^»,spin/-^orb (1/3 — i»,spin/iorb) 


a/aR 


CoRoT-1 


7.121 


31.657 


10.796 


0.612 


0.336 




1.667 


CoRoT-2 


2.208 


2.225 


3.850 


0.844 


0.186 




2.749 


CoRoT-3 


0.515 


0.650 


2.176 


1.272 


0.125 


0.208 


13.639 


CoRoT-4 


1.010 


1.010 


9.246 


0.998 


0.366 




6.141 


CoRoT-5 


14.930 


14.930 


59.986 


0.555 


0.113 




2.632 


CoRoT-6 


0.767 


0.767 


6.821 


1.205 


0.091 


0.242 


9.860 


CoRoT-7 


0.000 


0.000 


0.000 


0.149 


0.000 




2.764 


GJ1214 


3.385 


3.385 


5.336 


0.759 


0.698 




2.885 


GJ436 


17.070 


17.070 


45.131 


0.524 


0.130 




3.950 


HAT-P-1 


3.432 


3.439 


15.308 


0.734 


0.356 




3.294 


HAT-P-11 


5.108 


5.108 


25.289 


0.671 


0.543 




5.549 


HAT-P-12 


22.093 


22.093 


70.911 


0.545 


0.071 




2.517 


HAT-P-13 


9.328 


9.328 


27.208 


0.574 


0.278 




2.816 


HAT-P-2 


0.708 


0.708 


3.988 


0.995 


0.192 




10.659 


HAT-P-3 


29.067 


29.067 


84.263 


0.593 


0.034 




3.546 


HAT-P-4 


4.772 


4.772 


14.622 


0.712 


0.671 




2.722 


HAT-P-5 


8.142 


8.142 


22.702 


0.623 


0.150 




2.987 


HAT-P-6 


2.205 


2.205 


8.488 


0.848 


0.585 




3.506 


HAT-P-7 


11.113 


-11.124 


24.490 


0.553 


0.241 




2.804 


HAT-P-8 


2.004 


2.004 


6.949 


0.871 


0.602 




3.273 


HAT-P-9 


1.425 


1.425 


5.610 


1.036 


0.971 


-0.638 


3.055 


HD149026 


3.881 


3.968 


11.160 


0.868 


1.270 




4.094 


HD17156 


2.286 


2.294 


48.555 


0.945 


0.025 




20.702 


HD189733 


5.269 


5.270 


11.684 


0.722 


0.113 




2.809 


HD209458 


3.526 


3.537 


12.429 


0.731 


0.379 




2.799 


HD80606 


0.392 


0.651 


43.714 


1.054 


0.011 


0.323 


71.577 


Kepler-4 


10.631 


10.631 


34.186 


0.822 


2.159 




4.837 


Kepler-5 


5.325 


5.325 


18.893 


0.671 


0.208 




3.889 


Kepler-6 


7.236 


7.236 


23.451 


0.610 


0.315 




2.698 


Kepler-7 


4.543 


4.543 


22.194 


0.746 


0.816 




2.746 


Kepler-8 


2.034 


2.281 


7.158 


1.021 


1.273 


-0.939 


2.567 


OGLE-TR-10 


6.380 


6.380 


19.725 


0.631 


0.285 




2.698 


OGLE-TR-111 


0.000 


0.000 


0.000 


0.632 


0.000 




3.670 


OGLE-TR-113 


0.000 


0.000 


0.000 


0.575 


0.000 




2.340 


OGLE-TR-132 


0.000 


0.000 


0.000 


0.439 


0.000 




2.328 


OGLE-TR-56 


18.964 


18.964 


22.979 


0.489 


0.207 




1.734 


TrES-1 


10.363 


11.966 


31.397 


0.670 


0.065 




3.325 


TrES-2 


20.531 


20.787 


50.730 


0.613 


0.043 




2.957 


TrES-3 


20.960 


20.960 


27.380 


0.622 


0.039 




2.117 


TrES-4 


3.041 


3.060 


10.806 


0.834 


0.862 




2.410 


WASP-1 


5.259 


5.259 


13.252 


0.676 


0.538 




2.215 


WASP-10 


2.120 


2.120 


6.533 


0.988 


0.068 




4.431 


WASP- 11 


21.978 


21.978 


79.914 


0.631 


0.035 




3.449 


WASP-12 


33.152 


33.152 


36.094 


0.429 


0.185 




1.236 


WASP-13 


3.178 


3.178 


13.832 


0.779 


0.619 




3.169 


WASP-14 


5.964 


7.119 


13.481 


0.807 


0.050 




4.877 


WASP- 15 


4.984 


4.984 


18.676 


0.683 


0.520 




2.566 


WASP-16 


5.113 


5.113 


15.949 


0.695 


0.161 




3.746 


WASP- 17 


2.075 


-2.474 


7.755 


0.999 


1.228 


-0.894 


2.033 


WASP-18 


5.958 


5.958 


5.591 


0.713 


0.100 




3.522 


WASP-19 


14.951 


14.951 


11.759 


0.513 


0.244 




1.296 


WASP-2 


0.000 


0.000 


0.000 


0.578 


0.000 




2.815 


WASP-3 


2.673 


2.767 


4.945 


0.812 


0.612 




2.589 


WASP-4 


16.469 


16.469 


22.077 


0.566 


0.087 




1.852 


WASP-5 


9.151 


9.151 


14.870 


0.614 


0.134 




2.760 


WASP-6 


9.348 


9.348 


31.431 


0.629 


0.109 




2.717 


WASP-7 


0.742 


0.742 


3.677 


1.221 


0.978 


-0.644 


5.841 


XO-1 


10.799 


10.799 


42.559 


0.696 


0.051 




3.747 


XO-2 


13.409 


13.409 


35.080 


0.566 


0.121 




2.977 


XO-3 


1.827 


2.297 


5.811 


0.872 


0.166 




4.903 


XO-4 


2.159 


2.159 


8.966 


0.826 


0.382 




4.306 


XO-5 


18.603 


18.603 


78.035 


0.680 


0.030 




4.455 



